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ABSTRACT 


Timing and power dissipation are the two important perfor- 
mance indices of an MOS digital circuit. With this in view a 
program to simulate I'S®S/LSI circuit, named MOSIMR, has been 
developed. It provides the user with initial logic states and 
timing waveforms of requested node voltages and element powers 
(as well as total power dissipation) over the simulation period. 

The main focus is on minimizing CPU time by exploiting latency 
of large circuits. The given circuit is decomposed into indepen- 
dent subsystems by node partitioning. A linked data structure of 
the circuit topology is constructed to order the processing of 
node partitions in accordance with signal flow. This way a sele- 
ctive trace can be conducted through the circuit to take full 
advantage of circuit inactivity. 

Any change in input travels as a wave in a large circuit. The 
trace algorithm is capable of dynamically detecting this signal 
wave, so that only those nodes which are encompassed by the wave 
profile are analysed on a signal flow path. This results in 
further exploitation of latency. 

idth the above scheme of network - ordering, nonlinear Gauss 
Seidel (GS) relaxation method is employed to evaluate the active 
nodes only, in each time step. One iteration of relaxation, 
containing one iteration of Newton Raphson, is performed per 
simulation step. 

Results checked against SPICE 2G show reasonable accuracy and 
remarkable reduction in analysis time with increasing circuit size. 
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CHAPTER 1 
IMrRCOUCTIC»l 


Circuit analysis as a laeans to verification Of a design 
got a new dimension after tte advent of integrated liircuit (IC) 
and its rapid evolution over the past two decades. The tra- 
ditional papex-and-pencil analysis is a sound methtwji but it 
imposes severe lijaitations on the size and type of circuits 
that can be analysed economically. Large linear (say those 
containing > 50 elements) [b] or even small nonlineat circuits 
axe seldom analysed exactly. Instead engineers often rely 
heavily on intuition to obtain approximate solutions. Inva- 
riably the final verification of circuit performance ts obtained 
by breadboarding, where there is also scope, though limited, 
to optimize circuit performance by a trial and error adjustment 
of element values. 

But all these prove grossly inadequate for the analysis, 
design verification and optimization of integrated circuits, for 
it is impossible to duplicate the IC with discrete components. 
Certainly parasitic— effects and element-hatching characteristics 
between integrated devices are not reproducible in discrete form 
besides tolerance analysis etc, can not be performed with 
full freedom as we can scarcely control the device parameters. 
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A computer program to simulate the circuit can on the 
other hand build up circuit models of integrated devices to 
any degree of accuracy and provide control over all the device 
parameters. In fact it has proven to be the most effective 
tool for integrated circuit analysis, by which much more 
accurate results are obtained at a cost which is often a small 
fraction of the bread boardii^ cost [5], 

’with increasing compaction large scale (LSI) and very 
large scale (VLSI) integrated circuits are being designed and 
fabricated, A high premium is put on ensuring the correctness 
of design of these complex circuits before fabrication, 
because even a minor modification at later stage will result 
in repeating the expensive and time consiming operations of 
mask making. 

Concerted research efforts have been organised for the 
last two decades by different groups of people resulting in 
several general and special purpose circuit simulation programs 
some of which are enumerated below. 

1,1 GiJCRAL PURPOSE CIRCUIT SIMULATORS: 

Simulation program with integrated circuit emphasis 
(SPICE) is a well known program of this class. Originally 
developed by Dr, Lawrence Nagel, SPIOE was later extensively 
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modified in version SPICE IL by E. Cohen and D,0. Pederson and 
finally cot augmented to version SPICE-2G in 1981 by A.R. 

Nsivton, D.O. Pederson and A. Sangiovanni-Vincentelli in the 
Electrical Engineering and Computer science departnjent of 
University of California, Berkeley [8] • 

SPICE is a full fledged circuit simulator using the 
conventional method of analysis and contains circuit models 
for semiconductor devices like diode, BJT, FET and MQSFET for 
application in a wide variety of circuits, 

1.1.1 Anatomy of a Conventional Circuit Simulator: 

Most gene ral-fxir pose computer simulation programs have 
five main stages [5j which are elaborated below: 

1. 1.1.1 Input Stage : 

Here network topology and element specifications are 
fed in as input data to the program. Preferably this input 
should be in a format free user language which can be translated 
by a language compiler inside the simulator, into simulation 
tables. 

Error detection features (e,g, detection of loop of 
voltage sources, cutsets of current sources) are incorporated 
v/ith helpful diagnostic messages to make the program user 
oriented. 
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...ultiple level nesting of circuits, through subsciicuit 
definition by the user, is allowed to facilitate concise 
description of repeatitive structures. 

For a good user friendly program coding for this first 
stage turns out to be nearly 40/- of the whole program. 

1,1. 1.2 Device model retrieval and replacement stage ; 

The simulator has a built-in library of device models 
(e.g. diode, BJT, FET, MCSFET). New models defined by the 
user, as subcircuits, are also stored efficiently. 

The user has control over all the device paran^ters in 
the model. Any parameter unspecified by the user is assigned 
a default value by the program. 

An efficient retrieval scheme is kept for quick repla- 
cement of a device, when it appears within a circuit, by its 
equivalent model. 

The first two stages simulate the circuit in question 
and get it ready for analysis. 

A general purpose simulator can do three kinds of 
analysis i) DC-bias point analysis on linear and and non-linear 
circuits, ii) AC small signal analysis on linear circuits, 
iii) Transient analysis on linear and nonlinear circuits. 
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1 *1.1.3 Equa ti on formulation stage ; 


kind of analysis to be done 
Kin d of Ana lysi s 

i) Linear AC (frequency 
domain analysis) 

ii) Linear transient 
( time domain 
analysis) 

iii) Nonlinear d.c. 

(resistive analysis) 

iv) Nonlinear transient 
(dynamic analysis) 


summarised below - 

Equations formulated 

Linear nodal equations 
linear hybrid equations. 

Linear state equations 


Nonlinear nodal equations, 
nonlinear hybrid equations. 

Discretised nonlinear nodal 
equations 

Discretised nonlinear 
hybrid equations 
Nonlinear state equations 


Here circuit equations are framed according to. the 

as 


The equations generally in the algebraic-differential equations 
form represent the dynamical system either in the frequency or 
time domain. Prior to framir^ the equation for any kind of 
analysis however a EX; operating point calculation is performed 
to sat the initial condition* 

1. 1.1.4 Numerical solution stage ; 

The equations framed in previous stage are solved here 
employing different numerical methods depending on the types of 
equations, as tabulated below — 
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Typ e of equations 

i) Linear nodai/hybrid 
equations 


Method of solution 
employed 

Gaussian elimination, 

LU - decoraposition 


AT 

ii) Linear state equations By numerical evaluation of e 

(state transition matrix), 
Numerical integration (e.g. 
Backward Euler (BE), Tropezoidal 
rule (TR), Gear's variable 
order method (GE)) 


iii) Nonlinear nodal/hybrid 
equations 


iv) Nonlinear state equations 


Newton Raphson (NR) method. 
Piecewise linear method (for 
getting driving point and 
transfer characteristic curves) 

Explicit/implicit numerical 
integration 


1.1. 1.5 CXitput stage ; 


This final stage gives the results of analysis in any 
of the forms, summarised below, as requested by the user, 
depending on the kind of analysis. 


Kind of analysis 


i) Linear AC( frequency 
domain analysis) 

ii) Linear transient 

(time domain analysis) 


Forms of result obtained 
at output stage 

Frequency response »distortion 
analysis 

i>tep response, pulse response, 
periodic input response 
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iii) ‘..'onlineax DC 

(resistive analysis) 


Operating point calculation (DC 
bias). Driving point and transfer 
characteristic plots. 


iv) .lonlinear transient Transient response, periodic 

(oynarnical analysis) response for different v/aveforms, 

1.1,2 Disadvantages faced in Conventional Simulators: 


The conventional circuit sinwlators are versatile and 
accurate but they are cost effective for circuits containing 
only few hundred elements or less. For example SPICE IL can 
analyse circuits having at most 100 elements. In the quest for 
versatility SPICE incurs tremendous overhead (sec. A,l) in 
the first two stages (sec, 1.1.1). 


Sparse matrix method of solving the simultaneous eqns, in 
SPICE exploits sparseness of large circuits to save analysis 
time to some extent. But it is unable to detect and exploit 
latency (A, 3. 5) of large circuits. Besides SPICE operates in 
circuit level only, whereas in many cases it is possible to 
make with less accuracy than that provided by a full blovm 
circuit simulator. On the whole a conventional circuit 
simulator has proved to be expensive in terms of GPU time for 
the analysis of large circuits and has led to the development 
of special purpose simulators. 
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1.2 SPECIAL PLHPCbE SlkULATOftS: 

The special Purpose simulators axe meant for large 
circuits (e.g. LSI and VLSI) with dedicated use to avoid 
unnecessary overheads (Sec, A.l). The main aim in the design 
of such simulators is to improve upon CPU time. 

Simulators such as MOTIS [lO], SPLICE, DIAI4A, MCRO [llj, 
all special purpose simulators, in their quest for speed have 
rejected one or more of the principal features of conventional 
simulators. 

Fig, 1.1 shows the relative positions of these simulators 
on a speed/accuracy scale [l5J. On one extreme are the 
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1.2.1 Logic Simulation: 

A foremost example of logic simulator is Logic Analyser 
:or i.iainrenance Planning (LAAiP) [9], a digital logic-simulation 
program developed at Bell Telephone Laboratories around 1970 
'With elaborate test generation capability for fault diagnosis, 

LAMP simulates the circuit at gate level (INJOT, Nnll), 

NOR etc.) functional level (register, memory cell etc.) or 
at hybrid level. This sort of macromodeling helps to speed up 
the Simulator (Sec. A.2), 

lAMP also applies controllability and observability 
relation among functional nodes (compare with Sec, A, 3,9) to 
yield a directed graph representation of the digital circuit. 

Any loop within the graph due to feed backs in the circuit are 
•broken' by addition of control and observation at appropriate 
places to remove conflict (compare with Sec, A, 3, 10) anc build up 
a link graph (Sec. /%,3,9.5). 

A link graph gives an ordering of the nodes (Sec, A. 3, 9. 3) 
where a 'higher order’ node can be computed independently of 
any 'lower order’ node(s). 

This way LAi'jiP implements a very important concept of 
partitioning and ordering a network wrtiich helps xo bring dovm 
'-'PU time remarkably (Sec. 1.2#2,2).. 
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.^lso;the final link graph structure establishes a 
precedence order of signal flow and it becomes possible to 
conduct a selective trace along the signal flow path and 
exploit lateney (^aec. A. 3,5) of large circuits, with all 
these new concepts brought in LAMP is capable to simulate 
circuits containing 520CX) gates, economically on a IBM 370 
model 168 machine [9], 

SALOGS IV is a logic-timing simulator which models 
logic-gate inertial effects as pure delays rather than 
producing merely the docic levels. But it does not provide 
the analog waveforms. It is some ten thousand times faster 
chan a circuit simulator and is capable of analysing systems 
of many thousands of gates [l5]. MA^OT also models inertia 
as pure delays but gives an analog slope for rising and 
falling edges [l5j, 

1.2,2 Timing Simulation: 

Besides logic simulation timing analysis has been found 
to be of great interest in digital circuits. It is here that 
the user can get a more accurate picture of the profiles of 
different circuit variables (e.g. node voltages, device 
currents and power dissipation etc.) against time, rather 
than the simplistic picture of logic levels plus discrete 
propagation delays provided by the logic simulators. 
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TL'ning analysis for large circuits is very expensive 
//ith conventional circuit simulators. Research efforts in this 
direction have yielded special purpose simulators which 
incorporate many time saving features of logic simulators like 
i) Alinimising overheads ii) Macromodeling, iii) Network 
partitioning and ordering iv) Exploitation of circuit latency 
etc. This way these special purpose simulators have been 
able to achieve on an average one order of magnitude increase 
in speed over conventional simulators, while maintainif^ 
accuracy of results comparable to that of standard circuit 
simulators , 

Thus timing simulation stands mid way between logic 
simulation and standard circuit analysis and proves to be the 
best compromise between speed and accuracy of analysis, 

1 . 2 .2 , 1 T iming simulation technique s : 

There are three kinds of timing simulators depending 
on the techniques they employ. All the three kinds use 
relaxation method,but differently and at different stages of 
the algorithm for transient analysis (Fig. 1.2), The general 
relaxation method of analysing a system of algebraic-different- 
ial equations have been elaborated in Sec. A.3 of Appendix A. 
Representative examples of these three classes of simulators 
for time domain analysis are 1) MOTIS - a timing simulator 
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2) SPLICE - a mixed mode simulator performing iterated timing 
analysis (ITA) and 3) RELAX - using waveform relaxation method 
to perform timing analysis. 

Timing simulation done in MOT IS as well as Iterated 
timing simulation done in SPLICE apply Gauss-Seidel (GS) method, 
(:3ec. A. 3. 2) in the nonlinear algebraic equation stage of 
circuit analysis, thus carrying out a (^uss-Seidel-^ewton 
(i.e. Newton Raphson (!fl) iterations included within each 
G-5 iteration) method to decouple and solve the circuit 
equations obtained by discretisation. In MOTIS only 

single iteration of G-S and f'lR are performed in each time 
step (Sec, A.3,3) whereas SPLICE 1,6 iterates these to 
convergence [l2]. 

Timing errors (refer Sec, 4,1.3 and Sec. A, 3, 4) which 
are incurred in timing simulators for non-unidirectidnal 
circuits can be avoided in ITA. But large number of strong 
feedbacks often slow down the ITA process considerably as 
too many iterations have to be performed before convergence 
could be reached. 

Unlike MOTIS and SPLICE, RELAX applies thb relaxation 
technique right in the differential equation stage of 
analysis. Naturally the unknowns here are waveforms (elements 
of a function space) rather than real variables. The trial 
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solutr.on to start with is the guessed waveform of the unknown 
ov^r tha entire period of integration 0 to T, Accurate 
comparison between waveform-relaxation («ia) method, as it is 
called and ITA have not been carried out yet [ilj, [l2j. 

3ut aoth have proven to be very accurate and having guaranteed 
convergence properties whan applied to the analysis of 
circuits, 

1,2, 2, 2 Comparison with general circuit similator ; 

All the three classes of special purpose timing simu- 
lators out perform the general purpose circuit simulators 
speed wise, mainly due to the following reasons: 

1. Circuit deccxnposition by assignment partitioning (effected 
by the application of relaxation method) greatly reduces the 
complexity of analysis is (Sec, A, 3,1). A circuit which 
is not partitioned has complexity of analysis proportional 
to n where n is the size (i.e, no, of unknown variables) 
of the circuit and ai is a niraber ' depending on the sparsity of 
the circuit, o is small (large) if the number of inter- 
actions among unknown variables is small (large), m is 
usually between 1,2 and 2 [llj. 

But with decomposition, the complexity of analysis 
is usually proportional to n. Thus circuit decomposition 
is a major factor in reducing analysis time, specially in 
large circuits. 
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2, Circuit partitioning breaks up the whole circuit into 

subsystem blocks which can be analysed independently* A 

scheduler built on the same line as in logic simulators, 

order the computation that helps to speed up convergence 

it 

if ‘S-S relaxation is in use (See. A. 3, 4), As Vi-ell as/builds 
up a link graph structure of the subsystems and allows 
most naturally to schedule a method to exploit latency 
and avoid unnecessary computation of subsystems which axe 
inactive in a particular time step. As a matter of fact 
relaxation methods are ideally suited to exploit the latency 
(Sec. A. 3.5) of the circuit under analysis. 

3, Finally the special purpose programs built for specialised 

applications can avoid overheads (Sec, A.l) incurred by 

The speed plus accuracy 

general purpose simulator ./achieved by timing simulators is 
possible because they punch the scheduling techniques of 
logic simulators for ordering the computations of subsystem 
blocks, vfith the full fledged circuit analysis of conven- 
tional simulators applied to the computation of individual 
subsystem. 

1.3 TBih DOMAIN AimYSiS OF MC^ LSI CIRCUITS; 

The computation cost for tijae-domain simulation of 
large circuits is prohibitive. The cost lies principally in 
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function evaluations for the Jacobian. Again the whole of the 
Jacobian (coefficient matrix) must be memory resident if 
equation solution efficiency is not to be impaired [15]. Despite 
falling cost of computing memory, this latter factor remains a 
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ccnv3ntional circuit simulator, FOfy»i, which includes time 

required for discretisation and linearisation, is found to 

time fox 

dominate over SOLVE, which is the/solution of the set of linear 
equations, only in case of small circuits (havin-Q < 50 nodes). 
Special purpose simulators have come up to face the problem 
of transient analysis of large nonlinear circuits which appears 
so in surmountable with standard circuit simulators, 

MC6IMR is such a special purpose simulator desigr«d 
with following aims and implemented in FORTRAN 10 on a DEC 1090 
system, 

1.4 AIMS IN DEVEL0PIM3 MOSIMR: 

Aim here is to develop a special purpose transient 
analysis program for digital J4AiO&/LSi circuits. Speed and 
power dissipation are the roost important factors considered for 
the evaluation of an mOS LSI design. MOSIMR is meant to 
provide timing waveforms of voltages at nodes of interest and 
power dissipation in elements of interest, with reasonable 
accuracy. 

Focus is on improving speed of analysis and minimising 
storage requirement, to make analysis of large circuits with 
MOSIMR economical. 

The program is to compute initial condition based on 
a full-fledged logic simulation of the given digital circuit. 
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Hence it can be made useful as a mixed mode simulator that 
can perform as a logic simulator at t:0 and as a timing Simu- 
la tor subsequently. 

In the quest for speed, new techniques of dynamic 
detection of the traveling signal-wave profile are to be 
envisaged in program implementation, for fuller exploitation 
of latency. This also makes the program a useful tool as a 
signal tracer in large circuits. 

The program is also meant to be interactive. The course 
of simulation will be controlled and guided by the uses a 
large extent. 

1.5 THESIS CXn:LIHE: 

The following chapters dealing with the design, imple- 
mentation and performance evaluation of the NMC& timing 
simulator, MOSIMR, are organised as follows: 

Chapter 2 deals with defining the circuit analysis 
problem and formulating the theoretical scheme for solution, as 
followed in MOSIMR, It also mentions the assumptions made to 
simplify the solution scheme. 

Chapter 3 gives the ctetails of program implementation in 
terms of flow charts and explains the function of each program 
module, A glossary of variables and arrays used in the program 
has been appended here to make the flow charts understandable. 
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Chapter 4 gives the results obtained by analysing a fevi 
typical circuit configurations and compares the analysis time 
snd accuracy of the results obtained by with that of the 

results obtained by the standard general purpose ciicuit 
simulator SPICE 2G. 

Chapter 5 draws conclusion by pointing out credits and 
limitations of the program MC&IK1R, its various possible appli- 
cations and finally suggesting scopes of modification for its 
betterment. 

Appendix A is there to clarify some important terms and 
concepts through definition, brief explanation/ derivation and 
illustrative examples. 

Appendix B gives a complete manual fox ^CSIMR, to help 
and guide the user in building input data files for the simu- 
lator. It also makes the user aware and cautious of all the 
idiosyncrasies of the program. 



CHAPTER 2 


THBCBHTICAL FRAAEV^CftK OF MOSIMH 


2,1 DcFI:VlTION OF THE PROBLEM: 

The digital circuit to be analysed is described in 

terms of transistors/gates/functional blocks like pass tran- 
sistor/logic gate/latch. Each of these constitutes an element 
(Sec, A. 2). Interconnection of these elements are the fun- 
ctional nodes (or simply nodes). 

An n-node circuit having initial state is a dynamical 
system characterised by a set of mixed- implicit algebraic - 
differential equations of the form 

F(y»yf u) = 0 (2.1a) 

v(o) = V (2,1b) 

where, 

y(t)€ ^ is the vector of unknown variables, namely the 
n node voltages. 

y(t)CR^ is the time derivative of y at time t 
u(t)6 is the vector of r input variables at time t 
Yq € is the vector of initial values of y(t) at 
time t. 

'FI 

?:F. X R X R — > R is a continuous function. 


ifferential nodal equations 
, v(o) S 



,ue$cent state solution based on logic simulation 



«E. method 


Gauss-Seidel-'f»jewton iteration 
( 1 iteration of G S performed) 


Solution of Vi from the ith equation, corresponding 
to; ith partitioned subsystem (1 iteration of NR 
performed) 


Gate-advancement according to the scheduling scheme 


Is the end of trace reached ? 


Advance simulation time by one step 


.g. 2.3: General scheme of solution followed in i-<4uSHviA 
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Thus the problerr. of tiTO domain analysis of li'.OS LSI circuits 
taken up by hlCSllou-i boils cbwn to the solution of the dynamic 
equation (2.1a) with the help of initial conditions of 
eqn. (2.1b), 

The nodal equations for a MOS circuit with n-nodes can 
be written as follows: 

C(v,u)v = f(v,u) (2.2a} 

with initial condition 

v(o) = (2. 2b) 

where v€ is the vector of all unknown noae voltages and v^, 
the given initial values of v. 

. X 

u£R is the vector of all inputs 
f: X J^-R^ is a continuous function each component of which 

represents the net sum of currents charging the capacitor at each 
node due to the conductors and the controlled current sources 
(UOS devices) connected to the node, 

C: X R^ — R^^*^ is a symmetric matrix in which 

is the sum of capacitances of all the capacitors connected 
to node i, 

-C. .(v,u) ; i ^ j is the total capacitance between nodes i and j. 
V'/e can xev/rite equations (2,2) as follows, 

C(v,u)v - f(v,u) = Q 
v( c ) = v^ 


(2,3a) 
(2 -3b) 
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whicn are of same foiro as eqns. (2.1) and describe the dynamic 
system of the Mv.Oo LSI circuit which MOSikR seeks to analyze, 

2.2 AooUAi? Tiers : 

i) The circuit is considered to be at steady state at 

t = 0. Hence the initial state is obtained by a dc solution 
of the circuit, prior to its time domain analysis. 

ii) During computation of v^(t), except the voltage at 
node i, all other nodes are considered to be constant. 

This assumption has vital implication in terms of 
floating capacitors and node partitioning. 



C-) CO 


Fig, 2,1: Fs.eplacement of a floatii^ capacitor by equivalent 
grounded capacitances. 

Let us consider the floating capacitor between nodes 
(1) and (2) in Fig. 2.1a. i is tte current, sourced by node 1 
and sinked by node 2, which flows through the capacitor c. 


Heact = C Aj (V(l)-v(2)) 


(2.4) 
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,<h'?n v;e solve for V(l) we assume V(2) to be constant (equal 
to VC say) . 

jo the current charging the node capacitor at node l is 


iCi = C (V(l)-VC) = c 

SiDilarly the node capacitor current at node 2 is 

iCj = -C 

Which is equivalent to having grounded capacitors of value C 
at nodes 1 and 2 with currents ic^^ and ic^ charging them 
as shown in Fig. 2.1b, 


d V(2) 

dt 


This simple reduction of floating capacitors into 
equivalent grounded caf^citors helps to reflect floating 
device capacitances (e.g. at the input/output nodes of 

gates as shown in Fig, 2.2 and free the circuit from floating 
capacitor feedbacks. 



-in " Ss ^Qd ' ^out ^%d 


Fig, 2,2: Reflecting floating device capacitances at gate 
input/output nodes. 
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In the way as shown in Fig. 2.2 all floating capacitances 
axe replaceo by equivalent grounded capacitances xo diaconalise 
ant decouple eqn, (2.3a) as foliowsj 



^1 

^2 ^'^2 


f .(v. f 



u) 


* 



= f (v ,d ,u) 

UB 


(2.5) 


Except Vj^(the node voltage to be solved from the ith eqn.) 
in eqns. (2,5) all other node voltages involved in the 
ith equation are considered constant (by assumption (ii)) and 
have been included in the vector d^. 

2.2.1 Validity of the Assumptions; 

Assumption ( i) remains valid for circuits which axe in 
steady state at t=sO which means currents through * on* pass 
transistors and node capacitor currents must all be zero 
ini daily. But any attempt to analyse circuit which is in a 
transient state will make assumption (i) invalid. 



.is sumption (ii) becomes clearly invalid if time step 
of simulation is too large. For its validity an upper bound is 
imposod on the simulation step size, 

2.3 Gl^EnAL SCIEJC OF SOUJTIONt 

The flov^ diagram in Fig. 2.3 sketches the general 
schen^ followed by MOSIMFi in solving eqn. (2.3), apparent 
from the flow diagram the scheme is similar to that followed 
by a noniterative timing simulator (Sec. 1.2. 2.1). Relaxation 
method employed is one of Gauss oeidal iteration at the 
discretised nonlinear equation level. 

But new ideas have been introduced in schedulii^ 
selective trace (Sec, A. 3. 7) based on signal flew to make a 
fuller exploitation of latency. 

The general scheme after it has reached the 3rd stage 
shovm in Fig, 2,3, consists of three major processes namely 
i) Discretisation ii) Assignment partitioning iii) relaxation 
process, 

2.3.1 Discretisation: 

Discretising the ith eqn, in eqns, (2,5) by Backward 
Euler method. 
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Y/aere i = 1 to n and h = time step of integration. 

2,3,2 Gauss-beidel-Nev/ton i4ethod; 

This process involves application of Gs.uss Geidel 
iteration at the non-linear equation level with an NR intexation 
nested within each GS iteration. 

2. 3.2,1 Assignment partitioning ; 

This is related to a system of simultaneous equations. 

Each unknov/n variable is assigned to an equation in which it 
is involved (Sec. A. 3,1), However no tv/o variables axe 
assigned to the same equation thus giving n disjoint independently 
Solvable equations. 


Physically this amounts to partitioning the circuit 
into n independent subsystems each having an unknown node to 
be evaluated. 


From graph theoretic point of view assignment parti- 
tioning will be recognised as breaking of the subgraphs 
to resolve conflicts (Sec, 2.5), 

Assignment partitioning of eqn, (2,6) is performed 
using nonlinear Gauss-Seidei algorithm as follows; 


t+1 


11 


f t-M 

>1 • 


t4*l 


U 




(2.7) 


where i =» 1 to n 





is ‘Che unknown voltage assigned to the ith partitioned 
su osysteni 


d S, 


Ui - 2^ 

where is the set of nodes to be evaluated in tipve step 
^i 

t-f*l, after node i, 

and is the set of nodes already evaluated in tfl, 

J. 

before node i. 


and S 2 depend on the scheduling algoritlir. that determines 
the actual sequence in which the equations are solved. 

Except v^ all other variables have been included in 
and d in eqn. (2.7). These are the decoupling inputs 

j-i 

of the ith subsystem (Sec. A, 3,1) and by virtue of asfumption 
(ii) in Sec. 2.2 these are considered constant, while is 

computed from 2,7 in the time step t+1. 

2. 3. 2. 2 Newton-Raphson (MR) method s 

In order to solve from eqn. (2.7) I'ft method is 

applied to linearise the nonlinear current function f^ as 
follows: 


From eqn. (2,7) we have 
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fence root of the eqn. 


'i = 0 


{ 2 . 9 } 


the required solution of in tin^ step t+1, Now 




u, * °1 • ^ f 


( 2 . 10 ) 


ii + 1 . ^ t ■'■ 1 

where F’(v^ ) is the derivative of F,w.r.to v^ 

t+l 

and f* is the derivative of f-w.r. to v, , 


Hence by H.R. foxTiiula 


= F,(v7" 


)° + 


F!(v ^■'■^)i ,, 


( 2 . 11 ) 


114*1 1 

where ) is the solution 

t+l o 

and (v^ ) is the trial value 

Taking the previous tin^ step value of v^^ as the trial value 
and the value of v^ in this time step, as the solutionj 


/ t-i-l\0 t _ _ j / t+lil 

(v,- / = V. and (v. ) = v. 


( 2 . 12 ) 


ence from eqn. (2.11), and (2,12) we have 


Fj^(v7^) = Fj^Cv.^) + (v7^ - Vj*) F^(v7^) 1 ^ 


( 2 . 13 ) 
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Let ^y. 
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T.-i-n using eqns. (2.8), (2.9) and (2.10) ec.n. (2.13) gives 


A V. 


fi(v^ , 


t+l 


t+1 


U; 


,U 
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) 


- -Uj; 




(2.14) 


where is a nonlinear function of several node voltages 
in general and denotes the net current flowing into the node 
capacixance f|^ is the dynamic slope of the node capacit- 

ance current vs, node voltage and hence denotes output condu- 
ctance at node i. It is equal to the algebraic sum of the 
conductances of all devices having output at node i. Equation 
(2,14) is used by MOSiliiR to compute from v^^ (refer 

Sec. 2.4), 

2.4 VALIDITY OF Ti-fc bOLUTIONi 


Fig. 2.4 shows a simple case where i is the node to be 
computed with 0^^^ as the corresponding node capacitance and 
, the capacitor charging current flowing in from the 
inverter output. 


For the solution in eqn, (2,14) obtained by a single 
iteration of M.R. to be asymptotically correct, f|^ the 
slope used in eqn, (2.14), must be redefined. 
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Operating point 



= I driver - I load 


Figure 2,4 

In the limit h-^o egn. (2.14) reduces to ^ 

which is merely a difference form of eqn. (2.2a). 

In the other limit ® eqn. (2.14) reduces to 


AV. = f./f] 

1 1^ i 


(2.15) 


Fig. 2.4o shows the V~I characteristics of the load 
and driver devices of the inverter in Fig. 2.4a. Approxima- ^ 
th- driver curve by a line joining the operating point 

P and the final point F we have f, linearly varying with V 

* » , * 0 

v/ith a slope 



*'.edefining f‘ as fiG ,in equation (2.15), v/e get 

“ ’^of as ■ 
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Thus with the new definition of f* asyiriptotic validity 
of --t.n, (2.14) can be ensured ^though large errors in int.ir- 
riec’iute poin-s can not be avoided. is not easily computable 

huMCO in /.iOTIo instead , the line joining the operating point 
P with the origin (having a static slope S in Fig. 2.4b) is 
used to approximate the driver curve. This incurs greater 
error in intermediate points and results in some asymptotic 
error also as it ends up with the final voltage (refer 
Fig, 2,4b) instead of 


Fig, 2,4b) of the driver curve. This gives less error in 
intermediate points than the method using static slope and is 
specially accurate for small values of h. Hov«ver, it has 
no asymptotic correctness and may lead to total instability 
for large h. Therefore for eqn, (2.14) to remain valid, an 
upper-bound over h has to be imposed. 

2.5 imE PARTITIONIIC: 

Node partitioning as done by Gauss Seidel method in 
i'.iCwLViR is based on controllability considerations (Sec, A. 3, 11) 


concerned node as its output node, 
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.'.Col.ViR does not allow wired connection of outputs and 
v.hsrever there are noxe than one controls fannin^; into- a 
a node it is due to presence of pass transistors (as at node 
h in .-ig. 2,5c) producing conflict (Sec, A*3.9,'*)» IXjSV.'ui 
resolves such conflicts by foinraing a fan-in table (A, 3, 11), 

A partitioned subsystem centred around such a conflicting 
node contains the concerned node and all the fan-ih elements 
(A. 3. 11,1) of that node, v#iich have to be evaluated in order 
to evaluate the node in a particular tin» step, 

2.5.1 Tackling maC Subgraphs of Different Kinds (Sec. A,3,9,3)i 

Figs, 2.5a - 2.5c show the typical cases where digital 
circuits form MSC subgraphs. In each case MOSI!*iri resolves the 



Fig, 2.5(c) 
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In Fig. 2.5a presence of the floating cai^citor 
causes M to depend on a Iwer order (Sec. A. 3. 9. 3) node Nl. 
5ut in *V.OSl}.ixi these two nodes are decoupled by xecioving tiie 
floating capacitor and replacing by equivalent grounded 
capacitors reflected at the input and output. This resolves 
the conflict in computing N. 


In Fig. 2,5b computation of N depends on the lowr 
order node Ml, due to a feed back. 3ut MOSIMR uses the voltage 
value of N1 of previous tirag step t to cccipute N at t+1 
(Sec. A.3.4), This introduces a minimal error of 1 

step, generally, w^en using G-S method. 

One limitation of MOSIMR will be pointed out in this 
context. As the scheduling algorithm ®aployed in MC3SlMi is 
unable to distinguish between a forward path and a feedback 
path, this method of usir^ the previous time step value of 
the controlling node, tfitfien it is still to be evaluated in the 
present time step, is applied blindly, in Fig. 2,5b fox 
example if inputs are scanned in the order I^, I3 and a 
change is monitored in I2 the scheduler conducts a trace 

through nodes 4,6,7 and 8. When ccme to evaluate 6 i.e, 
node !l at t+1 it is found to depend on nodes 4,8 and 5, Of 
the three controlling nodes both 8 and 5 are unevaluated at 
t+1, of these two 8 can not be evaluated unless 6 is 


comput-d, Ho%¥ever, 5 can be computed independent of 6 as 
it C03S not belong to a feedback path. It is therefore 
physically possible to compute node 5 at t+1 before 6 is 
t-valuated. But as the scheduler in is unable to tell 

betv^een a feedback path (such as node 8 in Fig. 2,5b) and a 
normal forward path (node 5 in Fig, 2.5b) both nodes are 
treated equally by taking their previous time step values to 
compute 6, thus incurring further timing errors. 

In Fig, 2,5c node U depends on the lower order node 
N1 if the bidirectional switch (the pass transistor) is *on*. 

The ensuing conflict is resolved by kOSlMR durir^ node 
partitioning process. The partitioned subsystem assigned to 
evaluate N includes the invertor as well as the pass transistor 
(with node Ml)as fan-in ele:.)t.ivis (sec, A. 3, 11), 

m 

2.5 i'J-ThORK CBDERIMl: 

Node computation ordering is based on observability 
relations of nodes (Sec. A. 3, 9), This relationship is built 
in the program MOSlAiE in form of fan-out table or linked 
slemc-nt list [l9j fran the circuit description (Sec, 3,2), A 
selective trace is conducted through the circuit along the 
signal flow paths mth the help of a scheduler. This v;ay 
inactive /latent parts of th& circuit are excluded from analysis. 
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Thf node computation sequence is determined by a 
circ’jiar quc-ue [l9] which acts as a scheduler (ooc. .i.,3,3) 
in MCoLVii, It can be visualised as service queue in which 
the* ’front' pointer directs towards the node partitionCs) 
to be serviced presently, ’Rear', is where the pointer to 
the subsystem affected by the computation of the current 
node is pushed in. Between the front and rear pointers the 
subsystems remain lined up in the queue to get serviced in 
the 'first come first served* basis. 

Need for efficient utilisation of memory space has 
prompted the use of a circular queue in place of a linear 
queue. As more and more nodes/ elements get serviced the 
front pointer in Fig, 2,6 move down, leaving behind a space 



ear 


Fig, 2,6: The circular scheduler queue 

that can be used (shown dotted in Fig, 2.6), So when the 
rear pointer reaches the fictitious end { end of the linear 
array simulating the queue) the queue is wrapped round to 
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stoxe further elements which line up to get serviced in the 
specs left behind by the progressing front pointer. If under 
this circumstance the rear pointer tends to overtake the 
front pointer the queue is recognised to be full and circuit 
simulation is terminated,.. 

The circuit subsystems are visited in the order 
directed by tl^ scheduler and evaluated. This is called a 
circuit trace. The trace comes to an end when there is no 
element left to be serviced in the queue (i.e, when the queue 
is empty with the front pointer tending to overtake rear 
pointer). 

2.7 SIMJLATION STEP: 

The entire time domain analysis is broken up into a 
number of simulation steps upon discretisation. In each step 
all the starting nodes are scanned and updated. If any such 
node is fcund varying it is followed up by a trace directed 
by the scheduler along the signal flow paths. Coming to the 
end of trace (indicated by Q empty) next starting node is 
taken up. i/hen all the starting nodes teve been monitored 
and corresponding traces performed, the analysis is advanced 
to the next time step, : 

The time step of simulation h, is kept as an user 
adjustable parameter, chiefly determined by the node capacitance 
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value. Due to the constraint put on h in Sec. 2,4 and Sec, 
2.2.1, h must be chosen s;:.all enough to. keec^ 3. ./h ratio 

Xi 

sufficiently large. The best value of tir:e step found fox 
;.;OSi;.IR and kept as default value is o.l nsec, vvr.en the node 
capacitances are of the usual oxter of 0,1 p?. 

The main aim in developing i.CoI;.Ti. is improvement of 
speed of analysis v/ith aminimai stcrage xequiremx-nt . However 
if accuracy is jeopardised in the process, one way to trade 
accuracy with speed etc, will be adjusting h to ^even 
smaller value, 

2.8 INITIALISATION: 

Solution scheme in MOSL'.IR being noniterative tJw 
analysis cannot start with a guessed initial state. Rather 
an accurate initial state calculation is done based on a 
preliminary logic simulation of the digital circuit. 

The circuit is assumed to be initially at steady 
state and the initial state computation involves evaluation 
of all internal and external node voltages, device currents 
and power dissipation under quo scent condition. 

2.9 DEVICE MCDELS USED IN blQSI.VE: 

AS pointed out in Sec, 2.3 in order to compute a node 
all the elements included in the subsystem (i.e, elements 
which fan-in at the concerned node) have to be evaluated for 
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the current tin^ step, MCSIMR uses equational Eodel for tne 
devices constituting the elements. This avoids the irsnans*; 
space requireojent of tabular models [lo] but tends to 
increase computation time. In MOSIl.ft hov/uver very simple 
three terminal model (Fig. 2.7) is used to repr^-sent an 
NMOS transistor- 



Qjs Is 

Fig. 2.7: The simple model used in MOSIMR to represent 
an MiAOS device. 
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dielectric constant of oxide 
thickness of oxide 
width of channel 
length of channel 
s 


during evaluation of an element, Back gate bias effect and 

channel length modulation effect have not bev-n modeled. These 

however do not affect the result too adversely, 

2.10 SIGNAL~WAVE PROFILE DETECTION IN MOSBIR: 

2.10.1 Signal flow Patterns in Large Circuits: 








in the common portion of signal flow path* 


Fig. 2*8e 
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2.10«2 Latency on a Signal Flow Path: 

Like majority of physical systems i\,OS/LSI circuits 
demonstrate a lov; pass characteristic by taking a finite time 
to respond to any perturbation. Perturbations in the form 
of signals (i.e, varyir^ voltage levels) travel with a speed 
dependent on propagation delay of devices in iViOS/LSl circuits. 
This amounts to saying that even a node that falls on the 
signal flow path, as found from the link-graph, (Sec, A. 3, 9. 5) 
may remain latent in a particular time step if the corresponding 
signal wave profile does not encompass that node. This may 
happen in either of the two cases - i) The ^ve has already 
crossed the node, else ii) the wave is still to arrive at 
the node. In any case the active nodes on the signal flow 
path are those encompassed by the wave profile. In addition 
to selective trace M05IMR is capable of finding in any time 
step the start and end points of the propagating signal-wave 
profile and hence can detect the active nodes on the signal 
floiv path which are to be computed. Thus within a salectcd 
path again UCSIUR analyses only those nodes which axe active 
and makes a fuller exploitation of latency, 

2.10,3 Start and End Point Detection of a Signal ^avo Profile: 

profiles 

Figs, 2,8a - 2.8e show the typical signal-wav^ where 
y% = change in voltage, 1 = length of the circuit on signal 
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flow path , evidently 1=0 is the point of application of 
input. The actual algorithm follov/ed fox dynamic detection 
of beginning and u-nd points of a wave profile can be found in 
boc, 3.0,3 and it involves a node constancy/ variation flag 
table IFL which at the end of a time step contains the 
information v/hether a node has remained constant 

, IFL(i‘'^) = 0 , or has changed, 

IrL(N) = 1 . 

A change in node voltage level is recognised only when 
where £ = tolerance (Fig, 2,8a) 

In addition to ticfc step h, € is another variable 
parameter in MOSIMR which can bo adjusted by the user to trade 
accuracy with speed, C = ,01 mV gives accuracy comparable to 
a full-fledged circuit siraulator like SPIi2:, 

The wave characteristic of signal flow will be 
prominent in very large circuits such as MOS/Lbl, MOSIf.iR 
which is designed to take advantage of this feature and make 
fuller exploitation of latency, thus becomes very suitable 
for the simulation of large I'-aiDS LSI's, 



CHAPTER 3 


PROjRAiM BIPLEi4ENTATION DETAIL^ 

in this chcpter each aspect of rr.acromodelin^, circuit 
description, subcixcuit expansion, linked list structuring, 
partitioning, ordering, logic and voltage level inl,tiaii- 
sations, timing analysis and outputting will be taken up, 
revealing hovj these have been implemented in the program 

The program, written in FCETR.iii 10 (the language 
supported by the Fortran compiler of C£C 10) is of length 
approximately 32(X) lines. It has been tested extensively on 
DEC 1090 with different ft.£)S digital circuit structures of 
small and medium size, 

3.1 skeletal STRUCTIBE TIB PROGRAM: 

The program is divided into 6 functional modules as 
follows: 

1. Circuit description 

2. Initialisation 

3. Input 

4. Analysis 

5. (Xitput 

6. Control 
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Tabic 3,1 gives tne na.p.es of r.air. arid accci" aar.yir..; S'jcprogxsr.* 
-nd their respective functions . In*- :ollo**ving sections tai:e 
up the program modules one in; on .f ano •-•latorate the lx ir.ole- 
mentotion ivith che help of ilon cnarts, 

3.2 CIRCUIT DESCRIPT ICN tiCuULD: 

The circuit description module is cniefiy concerned 
with reading-in the linked subeixeuit and linked main circuit 
descriptions (Fig, 3.2), expansion of subeixeuits apd their 
linking with the main circuit environment (Fig. 3.3 end 
Fig. 3.4) and node partitioning by fan-in/fan-out table 
construction (Fig. 3.5). 

3,2.1 Reading-in Circuit Descriptions: 

The circuit description sy the user is in terms of 
circuit macromodels, defined as 'elements' in MOSIMR, in a 
fixed format specified in user's manual (Appendix B). An 
element can be a logic gate (e.g. invertor, two- input RAi'© 
and NOa gates etc.), a functional block (e.g. a latch) or 
simply a pass transistor, Rew functionsi blocks can be defined 
by the user as subcircuits. 

In main circuit description an element is specified 
by its name, number, inout/ouiput nodes, the name and number 
of the next brother element at each inpu": none and the narr.e(s) 



anc nu:.;ber(s) of the successor eiementCs) at the output 
nod6'’s) (exact format given in Appendix 3), o.pecificationc 
for an elemerjt of type K. and number are stored in i'lth ro'.v 
of the table meant for K., as element record. Any raain 
circuit element having a link vdth a subcircuit at node N 
is marked correspondingly with an *L* in the specification . 
Name and number of such an ele.nent are kept stored in r ovv 
N of LINK table fox future reference during linking of 
subcircuit with main circuit. 

The subcircuit description is analogous to that of 
main circuit. Only hero all node end element raemberin^s 
are local, and global linkages with the main circuit at® 
missing. All these therefore need processing during Smbcircuit 
expansion before final ele.nent records can be obtained. Hence 
the circuit description module upon reading— in the subcircuit 
description, stores specifications of each elemion t, along 
with its name and number, in each row of work area (i.e. IhK.% 
Table), temporarily. The pointers to the zone in ^i.e, 

tha I'/jKAR address where the description starts and the iU\<Jih 
address where the description ends), where the description 
of a subcircuit is recorded, along with the coricsponring 
sujcircuit name, are stored in a row of the pointer- table 


IPTR 
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3.2.2 Subcircuit Expansion and Global Linking: 

yhile reading the input data i^en a subcircuit call 
is encountered, the IFTR table is searched with the subcircuit 
naruG as key, to obtain the pointers to the zone in I.iiSiil 
where the subcircuit description can be found. Then slc-iw^nt 
records are fatched from the zone, processed and finally 
stored in appropriate elenant tables, sequentially. 

The processing involves 

1) offsetting local element nos, with the help of 
(refer Table 3,2). 

2) offsetting the local internal node nos. with the 
help of Jv.-i.XNC0 (refer Table 3.2) 

3) Replacing dimmy external nodes of the subcircuit with 
corresponding actual nodes provided by the user, 

4) Setting up links with the main circuit environment 
with the help of Llf^ table, 

processing 1-3 are elaborated in Fig. 3.3, As the subcircuit 
is expanded and new elements and nodes add up to main 
circuit, the element count vector ,..EL and node count variable 
are updated. 

processing 4, which involves linking of the subcircuit 
with the main circuit, can be of three types. 
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The suocircuit f?letr.ent fetched from t&ncorary store in 
and the main circuit el©i.ent kept xecordC'C in rov; ;i of 
LIIIK table axe to bo linked at node K, The possicilitios are: 

i) The main circuit element becomes the successor 
(oucc.) of the fetched elem, 

ii) The fetched element becomes the succ. of the main 
circuit element. 

In this case any future linking at node U would be 
’.vith the presently fetched element. Hence the present entry 
in row N of LINK Table has to be replaced with the name and 
number of currently fetched subcircuit eleinent. 

iii) The fetched end main circuit element can be related 
as brothers at the com.aon input node H. 

The direction of next brother linking is decided as 
follows - 

a) liihen the signal flows from the main circuit to the 

subcircuit, the fetched element is appended as next « 

chain 

brother (ND) at the end of the t«^of the main circuit 
element in LINK table. 

b) Otherwise the main circuit element is appended as h3 
at the end of the I® chain of rhu fetcned element. 
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Since in this latter case the fetched elr-ntnt cscor-.es 
the leading element (Sec. .-».3.6) of the overall :3-chain 
at node il, any future linking at N v*'ili have to xvfer 
to it. Hence present entry in rov/ h of LI?tv table ic 
replaced by the currently fetched subcizeuit eiement. 

All the possibilities (i) to (iii) are envisaged in Fig. 3,4 
and record of each linked and processed subcircuit element is 
finally stored in appropriate element tables. 

The end result is a linked list structure stored in the 
element tables, v/ith each eleracnt record holding a pointer 
(corresponding to each input/output node) to its M3/succ. 
element, 

3,2,3 Node Partitioning: 

The node partitioning algorithm by fan-in table 
construction has been elaborated in Fig. 3,‘ja, 3.f;b, The 
wired logic connection v/hich involves joining of tv/o or more 
gate outputs is not supported by MOSIMI. In fact tying up of 
tv;o or more drain/source nodes is treated as a conflict 
(Figs. 3,1a, 3,1c, 3, Id, 3,le) and is inadrdssible except in 
cases shown in Figs,3,lc-3.1e involving pass transistors. 

Here node N has, as opposed to ordinary case (Fig. 3, lb), 
more than one element controlling (Sec, A. 3. 9) it. Such 
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multiple controllability over H is resolved by luildinc a 
partitioned subsystem that consists of all che ele^Tents that 
fan-ia at iM, (Sec.A,3.11,l) , To compute .-1 it is necessary 
to evaluate all the elements in the partition assigned to node 

N. 


In AiOSIlifl the v/oxk space v/hich fcecor^s free 

after the expansion of all subcircuits, is used to store the 
fan-in as well as fan-out table for each of the conflicting 
nodes. As the only allowable conflict in ftXSL.il is associated 
with -j/S node(s) of pass transistors the node partitioning 
algorithm in the circuit description module takes up each 
pass transistor and explores the source node *N* fot a possible 
conflict, IFAN (N,3) stores the number of prcdececsoxs, as 
obtained from circuit description, at node N(refer Table 3.2). 
If at the explored source node the number of predecessors is 
> 1 a conflict is found to exist. The pulled end {considered 
as the drain terminal) of the pass rransistox as shown in 
Fig, 3,1c is always in conflict. The multiple controls 
(Sec, A, 3, 9) over these conflicting nodes are then tabulated 
in the 6th and 7th column of I'.wsAft by storing the names and 
numbers of the elements fanning-in (Sec, A. 3, 11,1) at 
(the node of conflict). Along side, the fan-out ele .lents 
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Fig. 3.1a: The inadrrdssible Fig, 3.1b: 
wired connection of Ordinary case 
gate outputs. with no 

conflict. 


Fig. 3,J.c: Gate 
output connected 
to D/S terminal 
of pass trans. 




Fig, 3,le; Bank of pas* transis- 
tors used as niultiplexei. 


(Sec. A. 3, 11. 2) at i\’ are also stored in the 8th and 9th 
columns of vjork area. Finally the pointer to fan-in /fan-out 
table location in the v«3rk area for a particular node N is 
stored in LliiC (:i,3) for future reference. 
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Construction of the fan-in table achieves node 
partitioning in the given circuit. The- node plus its fan-in 
elsiren'cs form a partitioned subsystem, 

3.3 INITlALIbATia^ MCCiULE: 

Initialisation of the given circuit is perfox ,ied in 
two sweeps. The first sweep carries out a full fledged logic 
simulation to determine the logic states of all the circuit 
nodes. 

The second sweep takes up each element and computes 
its output node voltage(s) and internal node voltages and 
device currents and the quescent element-power dissipation 
based on the input logic states, 

3,3,1 Logic Simulation: 

The nodes are initiated v/ith a don’t know state ’2’ . 

The user is required to specify the logic state of each input 
and the internal bistable elements. Starting a circuit 
trace from each of these nodes of known logic states, elements 
are visited in the order directed by the linked list, and 
their output logic states, canputed. 

The first three columns of. IWKAR wrapped around 
form a circular queue to schedule the servicing of elements 
(Sec. 2.6) and conduct an initialisation trace. 


AS the 
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logic state of a node N is determined, the leading number 
(Sec, A. 3, 6) of the chain of elements observing that node 
(Sec, A, 3. 9), is entered at the rear of the service queue 
to carry on the trace to the next logic level. On the other 
hand if the gate output remains undetermined after servicing • 
or is found to be already determined, further trace is 
terminated by not making any entry to the service Q, 

Naturally to determine' the ouput logic state of the element 
currently being serviced (in that, the front of the queue or 
any of its next brothers) either one or many of its inputs 
should have a dominant logic state (i.e, 1 for NOR gates, 

0 for NAND gates), or else logic states of all the input nodes 
should be known. The initial logic simulation has been 
detailed in Figs, 3,6, 3.7, 

Finally some floating nodes may remain (such as output 
node of an off pass trans,) which are now visited and the user 
is asked to specify initial logic states for these, to his 
liking. As each floating node is specified, a traqe is 
conducted from that node to initialise the subsequent levels 
(if these axe undetermined). Fig. 3,8 gives the details of 
floating node initialisation scheme. 

3.3,2 Node Voltage and Device Current Initialisation: 

At the end of logic simulation the user is asked to 
specify a number of NfviOS device parameters, interactively 
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through the TTY, These include threshold voltage (V^), 
surface mobility of electrons (jo.), thickness of oxide (t ), 
dielectric constant of the oxide ( , Vv/L ratio of load 

and driver transistor etc. Any parameter that remains 
unspecified by the user is assigned a default value by the 
program. With the help of these parameters intrinsic trans- 
conductance p is calculated for the enhancement mode driver 
device, the depletion mode load device and the pass 
transistor. 


Finally based on the known input logic states the 
simple drain current equations (2,16) are solved to get device 
currents, internal and output node voltages, and powej; 
dissipation in the element under quiescent condition, Follov/ing 
assumptions are made during solution of the initial element 
state; 


i) Input logic state 

1 

3 

0 


Corresponding voltage 
level. 


^DD = 

^DD- 

< 


5V 

= 


4V 


ii) Circuit is initially at steady state, implying zero 
current flow through 'on' pass trans. and hode 
capacitances. 
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iii) The saturation current CLS (computed from the device 
parameters) flov/s through the load transistor under 
steady state when the gate is ’on*. The drain 
currents are zero v;hen the gate is 'off*. 

iv) An ’off‘ pass transistor has zero subthreshold 
leakage-current. 

The actual voltage level for logic low (assumed < Vj) 
has no influence on restoring logic elements. However in a 
steering logic element like a pass transistor, the dtain and 
source voltages track each other if the element is »on*. If 
the input voltage, in such a case, is not known during evalua~ 
tion of the pass transistor, the output voltage remains 
unevaluated. This may happen if the input logic is ’O' and 
its voltage is not yet computed. MOSIi/iR assumes a default 
0.5V for the output, in such cases. But this approximation 
causes some error in the voltage initialisation (Sep* 5*2)*, 

3.4 INPUT iviODULE; 

The inputs can be either externally defined (EXDEF) 
or internally generated (INGEN). EXDEF inputs can be any 
arbitrary waveforms, values for v/hich are tabulated for every 
time step over the entire simulation period. This table of 
externally defined inputs is provided by the user fis input 
data to the program. 
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An INGEN input on the other hand is a regular input 
such as i) a continuous 1 phase/2 phase non-overlapping clock 
waveform, ii) A rising/falling edge iii) A positive/negative 
pulse or simply a iv) constant logic low/high level. These 
are inputs useful in digital circuits and are internally 
generated by MOSIMR based on user given specifications. The 
specifications include i) Rise/fall time ii) High/low 
period of a pulse etc, and are accepted from the user in the 
input module as illustrated in Fig, 3,10. 

3.5 ANALYSIS MODULE: 

After computation of the initial state the linked and 
partitioned circuit is subjected to the specified input signals; 
and its time domain response is analysed. The analysis module 
which follows after the input module does this task by condu- . 
cting a selective trace along the signal path and cofliputing ; 
only the active node partitions, l 

I 

3,5.1 Selective Traces i 

( 

The same circular service queue implemented on I^VKAR 
and utilised by initialisation trace, becomes the hiainstay of 
analysis trace as well. An analysis trace is initiated by 
a change in the external input which perturb the system that | 
is at steady state to start with. In each simulation step 
the external inputs are scanned and updated. Any change in 
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any input at node N greater .than a preset limit ^ , is 
recognised as a signal which will affect the elements 
observing (Sec. A. 3, 9) N, These are the elements obtained 
as succ, (or its next brothers) at N or as an entry in 
the fan-out table for node N, The leading element (Sec, A. 3, 6) 
of these is entered at the rear of the scheduler queue. This 
way the partitioned subsystems are ordered on the basis of 
observation (sec, A. 3, 9) by the scheduler, as each node is 
computed. Figs, 3,11 - 3.13 show the scheduling algorithm 
in detail, 

3,5,2 Computation of a Node Partition; 

’Front* of the scheduler queue points at the element 
being currently serviced. 

Servicing an element actually involves finding the 
contribution by the element in the present time step to the 
voltage(s) at its output node(s), N. The contribution is 
computed by solving the eqn, (2.14). Using current values 
of the decoupling inputs (Sec, 2, 3,2,1) (recognise this as 
the Gauss Seidel method), each element in the partitioned 
subsystem is evaluated (which involves updating of its 
device currents, internal node voltages and power dissipation 
stored in the real element table) and its contribution to 
node N, found out. Accumulating the contributions by all 
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elements which have control (Sec* A. 3, 9) over N, the node is 
updated for the current time step. For a node of conflict 
(Sec, 3,2,3) the elements included in the subsystem, are listed 
against the node in fan-in table. Figs, 3,14 and 3,15 give 
the details of node voltage updating in a partitioned subsystem 
with individual element evaluation, 

3,5,3 Signal Vi/ave Profile Detection: 

This has been elaborated in Figs, 3,12, 3,13, In 
essence it is based upon the following points; 

i) NCON is a flag storing the infomation whether the 
currently computed node was constant in the previous 
time step, 

IFL (NCBE) is the constancy flag of the currently 
computed node in the present time step, 

IFL(N) is the constancy flag of the starting node N 
from where the analysis trace has originated, 

ii) Let VOUT be the evaluated value in this time step and 
VNCOE (NODE) the previous time step value of the voltage 
at node ’NODE’, Then, 

= ABS(V0UT - VNODE (NODE)) 
is the change in voltage at NODE in this time step. 

If ^ >C recognise it as a signal at NODE 
/\ < ■£. recognise it as no signal at NODE 
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Based on this simple premise the analysis trace detects the 
starting and finishing points of a traveling signal wave in 
each time step, Vihile conducting the analysis trace, upon 
evaluating each node on the signal flow path, is checked 
against € , a preset limit. 

If A > €. the flag IFL(N) is checked and if it shows 
that the starting point of the current analysis trace is 
constant, in the present time step, at once the current node 
'NODE’ is recognised as the new starting point of the signal 
wave. Further progress of current trace is stopped and the 
name and no, of the predecessor element on the trace path at 
•NODE* is stored in IV^KAR (NODE, 10) - IVfiCAR (NODE, 11) to mark 
the starting point of fresh traces in the current and subsequent 
time steps. This way the analysis of leading part of a 
signal flow path, which is found inactive, can be skipped 
(refer Fig. 2.8b), 

In a similar fashion end of the signal wave can be 
detected when ^ <C (with starting node 'N* of the trace not 
constant) and rest of the trace, beyond the signal wave 
profile, being redundant, is abandonned (refer Fig. 2,8a), 

This simple scheme of signal wave profile detection, 
incorporated in MOSIMR, is useful in exploiting latency of 
very large circuits, where only a fraction of nodes/el^ents 
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remain encompassed by the propagating signal wave profile, at 
a given instant of time, on a signal flow path. 

3.6 OUTPUT MODULES 

Before the start of timing analysis the time step and 
simulation period are specified by the user in nano-seconds. 

The time step DELTA (0 < DELTA ;< 1 n-sec) is the one used in 
analysis for discretisation (by B.E. integration), whereas the 
printing/plotting step is computed in the output module as 
shown in Fig. 3.17. It is kept roughly around 1 nsec, for all 
choices of DELTA. 

At each printing/plotting step, the current value of 
the requested node voltages and element powers, obtained from 
VNODE and real element tables (refer Table 3,2) are printed 
and/or plotted against current simulation tirr^. The user can 
ask for at most ten node voltages and ten element powers to 
be outputted this way. 

Total power, if requested for, is calculated and 
printed/plotted similarly. The total run time required by the 
job. is printed at the end of the output file. 

3.7 CONTROL MODULE; 

At the end of simulation period the control module 
interacts with the user, offering him the options to i) COiTTINlE 
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ii) STOP or to iii) PRESTART further simulation. Thus user 
is allowed to steer the future course of the simulator 
(Fig.. 3.17). In case of CONTINUE’ ing^ the user is asked to 
specify the new simulation period (he is also allowed to 
specify inputs (not all) for the new simulation period). The 
present state of the circuit is preserved, and further 
analysis is continued from this state, as control is returned 
to the beginning of analysis trace. 

In the case of RESTART' ing, the control is re-turned to 
the beginning of Initialisation module. Present state of the 
circuit is destroyed and analysis starts by reinitialising 
the circuit, performing a logic simulation. 



Table 3,1; 


Main program MOSli\/iR uses following subprograms to implement 
the 6 functional modules. 

Name of the subprogram Functions performed 


SUBl 

SUB2 

Circuit description 

SUB3 

SUB 4 

module. 

LIMIT 

ITRACE 

LEVAL 

Logic initialisation 

VIM IT 

voltage/current 

initialisation 

INSPEC 

IMGEN 

Input module 

ATRACE 

^EVAL 

SBRTN 

eleval 

Analysis module 

PLOT 

Output module 


N.B.; The control module is included within the main 
program itself . 
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iable 3,2? A glossary of arrays and variables used for 
structuring cata in 

Name of Diinsnsions Functions 

array/ 

variable 

III 1 head pointer. Useful in diagnostic messages 

to locate error in input data, 

lUiUul 100 X 18 Viork Area, i) Temporarily stores sub- 

circuit descriptions (before the subckts, 
are expanded) - Circuit Description module, 

ii) Used as a circular eueue v/hile condu- 
cting logic initialisation trace - 

Initialisation module, 

iii) Holds the spv.-cif ications for the Ith 
input in lUivU: (1,13)- I7.iUtPi.(l,17) 

- Input module, 

iv) a) IV.ix-tIi(l, 1)- I'riCtfi(1, 3) used as a 
circular queue to conduct Analysis 
trace, 

b) I 1 , O ) — I vii1CiHH( I , 9 ) ■" holds 

fan- in/fan-out table for appropriate 
nodes. 

c) I.CsAR (I,10) — I.<i'CiR( I , ll) stores 
the name and no. of the predec, elem, 
at node I for starting a trace from 
the internal starting point I, 



fu T-i 
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l\!c 



Dimensions 


Functions 


IPTFi 


iv"ibL 


!■ OD 


10 X 3 


10 


1 


iii) LI:'IK(I,3) holds the pointer to 

the I;.KAfi rov; from where onvBrds 
fan-in/fan-out table is kept stored 
for node I. - Circuit Description 

module. 

Pointer table, i) Pointer to the stored 
description of each subckt. 

IPTrL(l,l) - stores subckt. name 
IPTii(I,2) “ holes the I0K.tFi address from 
where the description begins, 
IP1TI(I,3) - holds the address where 

the subckt. description ends, 
IMo. of subekts. that can be described 
should be < no. of rows of IPTR 

i.e. < 10 ~ Circuit Description 
module. 

Element Count, 11.1(1) stores the no, of 
elaraents of type I, present in the circuit. 

Node Count, No, of global nodes present 
in the circuit. 

Predecessor table, stores the predecessor 
element type and no,, for each node , 

IFM(I,l) - stores elem. type of the 
leading predec. element 

contd. . , 


IFAN 


100x3 
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of 

array/ 

variable 

Ei i 

EXTNOD 

VNODE 

CAP 

INVER etc. 


Dimensions Functions 

IFAi'J(I,2) - Elem. no. of the leading 
predec. elern. 

IF,^iN(I,3) - stores the no. of predecessor 
elements present at node I, 

- Circuit Description module, 

9x3 External link. This contains in each row 

a duimuy node no, and the name and. no, of 
.the subckt. elem, which has external conne- 
ction v/ith the main ckt. at the dummy 
node . 

9 This stores the actual nodes which replace 

the duOTiiy nodes of the called subckt, 

100 Node voltage vector. ! 

VNODE(I) stores the voltage of node I in 
the current simulation step, 

100 Node capacitance vector: 

C/\P(l) - stores the equivalent (device + 
stray) capacitance at node Vi/ith respect 
to ground, 

. 

50x6 Element table. Each row of an element 

table contains an element record. For j 

example element record of the Jth inverter 
is as follows; i 


c ontd 



Dimensions 


Functions 


Name of 

array/ 

variable 


IvElN'V etc. 50 x 5 


IFL 100 


InV'cH(I, 1) - stores the input node no. 
INVER(Iy2) - store a the output node no. 

IiWlR.(I,3) - INVEFi(I,4) - store the name 
and no, of the next brother element at 
the input node, 

IiWcR(l,5) - INVEFv(I,6) - store the name 
and no, of the succ. element at the output 
node. 


Real element table. This contains the 
real records of elements. For example 
racord for the Ith inverter is stored as 
f ollov/si 


, l) - stores the load transistor 

current in the present simu- 
lation step. 

kEINV(I, 2) - stores the driver transistor 
current in the present simu- 
lation step, 

FIEINV(I,3) - stores the power dissipation 

by the inverter in the current 
simulation step. 


Node constancy flag, 

IFL(I) = 0, v/hen node I remains constant 
in the present time step, 

contd. . 
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Name of 

array/ Dimensions Functions 

variable 


iFL(l) = 1, when noce I varies in the 
present time step* 

4 Temporary storage vector. These temporary 

storages are used to pass successor or 
next brother elements into Initialisation/ 
Analysis trace routines during scheduling. 

i NSTuP - simulation step counter 

NPEK - stores the specified petiod of 
simulation. 

The following arrays' and. variables are used in MLVAL anti ELEiVAL 
routines Vt/hich are respectively involved in evaluating a node and 
evaluating an element. 


1 3 o r NB 


NSTAP and 
NPER 


NOUT 2 


MODOui 2 


VOUT 2 


iiOUT(l) and i'>iuUT(2) store the two possible 
output nodes of the element being computed 
in ELdVAL routine, 

NODOUT(l) and ijODajT(2) contain the two 
possible output nodes of an element on 
signal flow path. These are evaluated by 
the i^VAL routine. 

The temporary store of the new voltage at 
K’ODOUT(I), after the node has bean updated 
by the NAVAL routine, for the current 
time step. 


contd. . 
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Name of 

array/ 

variable 


I£FL 


CUMUiC 


DELG 


Dimensions Functions 


2 Evaluation flag, 

IEFL(I)=1 if NCXJT(I) is already 

evaluated in the present 
time step, 

lEFL(l) = 0 if NOUT (I) is sti.ll not 
evaluated in the present 
time step, 

2 Cumulative capacitor current, 

CUMULC(l) cumulatively stores 
Contribution of each fan-in element 
towards the node capacitor at out^jut 
node NOD(Xrr(l). Finally it hol,ds the net 
capacitor current charging the node 
NODOUT(I), when all fan-in elemants at 
NODCXJT(l) have been evaluated, 

2 Element contribution to capacitor current 

DELC(l) holds the contribution of the 
element, being evaluated by EIJEVAL, towards 
the capacitor charging current at node, 
NODOUT(I). 


contd... 
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Name of 

array/ Dimensions Functions 

variable 


CUMULG 2 


DELG 


2 


Cumulative output conductance. CUMULG(l) 
cumulatively stores the contribution 
of each evaluated fan-in element towards 
the output conductance at node 
NODOUT(I). 

Element contribution to the oytput 
conductance, 

DELG(I) holds the contribution of the 
element being evaluated by ELEVAL, 
towards the conductance at node 
NODOUT(I). 
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CHAPTER 4 


results 

Different small and medium sized MMOS circuits have 
been analysed and the results compiled in this chapter. 
Comparison with results obtained by SPIC2i 2G shows relative 
improvement in simulation time, with increasing circuit 
size. The inaccuracies incurred by MCSIMR have also been 
made apparent by checking against the SPICE results^ 

In the voltage/power plots the points are joined with 
pencil for the monitored waveforms. The peaks (or humps) 
are marked with P's (red ink) and dips (or depressions) are 
marked with D's (green ink), showing their respective timings. 
Voltage levels of interest (such as final value reached or 
the tip of a hump etc.) are demarcated with blue lines. 
Finally any range of time/period, of interest (such as 
rise/fall time of an edge, an oscillation period or a pulse 
width etc.) are marked with black lines. 

4,1 TV^O PHASE DYNAMIC SHIFT REGISTER CELL? 

This is a typical example of a dynamic circuit, driven 
by clock. With no provision for subthreshold conduction of 
'off* pass transistors, MOSB-IR can not figure out the true 
performance of such a dynamic circuit for slov/ clocking. So 


V'le apply here fast clock and rapidly changing input to check 
the peri'ormance over 50 nsec,, which is much less than the 
usual period of refreshing (^msec.). 

The voltage at node 3 follov/s the rising edge at 1 
with a lag, to reach the final value ^ 3.75V. This causes 
the inverter output to fall at node 4. The rise time of node 3 
is 13 nsec, and fall rime of node 4 is 10 nsec. When phase 2 
of tihe clock at node 5 rises, pass transistor 2 turns ’on*, 
tending to equalise the drain, source potentials by redistri- 
buting the charges on the two node capacitances. By initial 
state assignment the floating source node 6 is at logic low 
(i.e, 0.5V in MOSIIVB and 0 V in SPICE) while the drain node 4 
is at -~^,8V, Hence the driven end 4 is pulled down 
momentarily by the floating end (which is at a lower potential) 
when the clock input at node 5 rises, to produce a depression 
at D. 

Comparing Figs. 4.1c and 4.le the timings and voltage 
obtained by MOSIMR are found to be in agreement with that 
obtained by SPICE. The element power dissipation as calcu- 
lated by MOSIIvlR over the simulation period have been printed 
and plotted for IN^^RTER 1 in table 4.1 and Fig. 4. id. A 

sharp rise in power dissipation occurs between 4 nsec, and 

inverter output makes a 

10 nsec, as the/transition from high to low in this period* 
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For analysing this circuit (v/ith 6 trans.) ratio of 
time taken by SPICE s time taken by MOSIMR = 4.79 

4.2 NMD L AT CHS 

This is a typical example of a tightly coupled feed 
back circuit yielding a bistable element. 

Initially node 3 is assigned ‘low* and node 4 'high* 
(i.e. reset state), A low going pulse is applied at the input 
node 1 while node 2 is held constant 'high'. As the flip flop 
sets, outputs 3 and 4 are monitored and the waveforms roughly 
agree with that obtained by SPICE. The dip D occuring in 
V(3) at 31 nsec, and a corresponding peak P in V(4) (the peak 
is missing in SPICE plots because of voltage scale cramping) 
are caused by the rising transition of the pulse V(l). Fig, 

4, 2d shows the total power dissipated by the Latch plotted 
against time. The power peaks at P^ and P 2 occur when both the 
NAND gates conduct under transient condition. The quescent 
power dissipation of 0.276 mw computed by MOSIMR, exactly 
coincides with that computed by SPICE (refer listing 4.2). 

For analysing the Latch (a circuit with 6 trans.) the 
ratio of time taken by SPICE : time taken by MOSIMR =* 4.08. 

4.3 MOS RING OSCILLATOR: 

This is a typical example of a circuit with a global 
feedback to yield an analog o-scillator out of digital gates. 
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1 he oscillation is triggered by a falling edge at 
input 1 and the waveform at the feedback node, 4 is monitored. 

The peaks and troughs of the periodic waveform at 
node 4 are 4.4V and 0.78V respectively in both M03livlR and 
SPICE plots (Figs, 4,3c and 4.3e). Hovt/ever, the period of the 
waveform obtained by SPICE is 4,8 nsec, while that obtained 
by iviOSlklR is 5,7 nsec. The discrepancy is due to the 
timing error, unavoidable in such nonunidirectional oircuits 
(Sec, A, 3. 4), The timing error can accumulate over the 
period of simulation and result in considerable deviations as 
clearly shown in Figs. 4,3c and 4,3e, The fcOSB/lR plot 
completely misses the trough obtained by SPICE due to 
timing error of 1.5 nsec, cumulated over the simulation period 
of 15 nsec. 

The waveform V(4) is near sinusoidal but with a slightly 
slower rise and sharper fall, caused by slow pull up through 
a load transistor of higher resistance and faster pull down 
through the low. resistance driver transistor in inverter 2. 

Power plot in Fig. 4,3d exhibits a period of oscilla- 
tion <^1.8 nsec, which is around ^ of that of V(4), 

For analysing the Ring oscillator (a circuit with 
7 transistors) over 15 nsec., the ratio of time taken by 
SPICE i time taken by MOSIMR = 7,77, 
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4.4 ONE -BIT FULL ADDER ViITH PASS-TElANSiSTOR CARRY CHAIN: 

This exarnplifies a typical combinatorial circuit 
whose response to the input v/aveforms, shown in Fig. 4.3a 
has been studied by monitoring the SUM (i.e. node 13 in 
MOSlivR circuit) and CARRY (i.e. node 9) outputs. The initial 
voltage assigned to V(9) which is at logic low is qV in SPICE 
and 0,5V in MOSIMR, Discrepancy K observed in V(9) (Fig, 

4,4d) is due to this difference in initial value of V(9). 
Further deviations can be noticed in the rise and fall slopes 
of the SUM and CARRY waveforms. However, the gross features, 
such as width of the SUM pulse (25 nsec, in MOSBfi and 27 nsec, 
in SPICE) and final voltages attained by SUM and CARRY in 
logic high/lov/ states, appreciably agree. 

The ringing observed at K in V(13) (Fig, 4.4c) is due 
to numerical over shoot (due to the sharp transition of Uie 
falling edge of the SUM pulse Sec, (2,4)), avoidable by 
lowering the time step of simulation. 

The rise in CARRY waveform starts earlier than the fall 
in SUM waveform due to propagation delay in inverter 3 
(Fig, 4.4a) . 

For analysing the Adder (a circuit with 21 transistors), 
the ratio of time taken by SPICE : time taken by MOSlMiR = 9.69. 
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Listing of the Input files for this circuit accepted 
by MOSB/iR has been given in Listing 4.4 to illustrate the 
user description of a circuit that does not contain any 
subcircuit. 

4.5 TWO-BIT JOHNSON RING COUNTER; 

This circuit has been chosen as a typical exa»nple of a 
finite state machine structure, it is a switched tail riiig 
counter with 4 states. Analysis starts from the reset state 
(the counter is reset by making the reset input R = ’!') 
recognised as Count = 0 and represented by the output states 
V(9) = '1' ; V(10) = 'O'; V(ll) = 'O' ; V(12) = 'O'. There- 
after the reset is removed by making R = 'O' and one cycle of 
R^-clock is applied at nodes 1 and 2 (with 0^ initially high) 
to effect transition to the next state i.e. Count = 1, 
(represented by V(9) = 'O' ; V(10) = V(ll) = 'O' % 

V(12) = 'O') while V(9) and V(lO) are monitored. 

In the plots obtained from MOSLViR and SPICE (Fig. 4.5c 
and 4.5e respectively) V(lO) is found to start changing 
earlier than V(9) as the count crosses over from '0» to *1*. 

This is due to the propagation delay introduced by NOR gate 2 
in the subcircuit-, BLOCK, between nodes 5 and 4. The rise/ 
fall times of the monitored waveforms (-»-'5 - 6 ns) and 
voltage levels (^0.58V when logic low and 5 V when logic high), 
obtained by MOSBvlE and SPICE agree very satisfactorily. 
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A hump P in V(9) occurs between 49 ns — 53 ns, synchro- 
nised with the rising edge of 02* It is due to the momentary 
lowering of V(4) by the redistribution of charge between the 
subckt. nodes 10 and 4 as pass transistor 3 is turned ’on’, in 
the subcircuit BLOCK. The voltage level and timing of the 
hump, given by MOSlklR roughly matches with that of SpICE. 

Total power dissipation in the circuit is plotted 
against time by MOSIMR in Fig. 4.5d. 

The input data for the circuit (which involves subckts.), 
as accepted by MOSiMR and SPICE are given side by side in 
listings 4,5a and 4,5b to make relative merits and demerits 
apparent. The forma t- less inputting of SPICE is more user 
oriented than the formatted input of MOSIftE. But absence of 
macro blocks, such as logic gates, in SPICE necessitates these 
to be described as subcircuits before being used in the main 
ckt. In some cases therefore, a SPICE input file may turn out 
to be more lengthy than that of MOSIk'iR. In the given example 
the input files are of roughly equal lengths. 

For the analysis of this medium sized circuit with 
52 transistors, the ratio of time required by SPICE : time 
required by MOSIMR = 13.6. , 


97 


4.6 CIRCULATING REFRESH REGISTER WITH R/W CONTROL: 

The refresh register is tested in the ’write’ mode with 
cototrol input (i.e. input 1) held at '1'. Floating node 3 
in subckt. CELL is assigned logic low' initially. Data input 
to the upper bank of registers, D1 (i.e. input 2) is held 
constant at 'O' and. the data input D2 (i.e. input 3j is held 
constant at so that only the bits in the bR cells belonging 

to the upper bank cnange with clocking, from 1 to 0» w/hile 
the lower bank remains unchanged. One cycle of the % phase 
clock is applied at nodes 12 and 13 while nodes 5 and 9 are 
monitored, V(5) changes from logic high to logic low and V(9) 
remains constant, at logic high, during this clocking, as 
expected. The gross nature of the waveform V(5) (nan^ly, the 
fall time, timing of the hump P or the steady value attained 
by V(5) at logic low ( /^ 0,7V)), given by MOSlklR tally with 
that obtained by SPICE, 

However, the circuits analysed by SPICE and MCSLVfi in 
this case though logically equivalent are not physically 
identical (refer Figs. 4.6a and 4.&b). In the circuit used 
by SPICE, the SR cell consists of pass transistors and 
inverters, whereas the SR cell, in the circuit used by MOSIMR, 
consists of pass transistors and inverting buffers . This 
results in deviations in the details of the voltage waveforms 
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V(5) (“‘Q* SPICE, using inverters in place of inverting 
buffers in the SR-ceils, gives a smaller hump at P compared 
to MOSIMR) . 

Though SPICE analyses a simpler circuit with 58 transistors 
and MOSIMR, a more complex circuit with 82 transistors, over a 
simulation period of 50 nsec., the ratio of time taken by 
SPICE : time taken by MOSlffi is 24,12, 

•* 

This large gain in the speed of analysis by MOSliOR is 
due to detection of latency of the lower bank of regitters, 

MOSIMR constrains its analysis to the upper bank of registers 
which undergo change and proceeds along the SR chain only 
up to that point which is active. Thus leaving out a major 
fraction of the circuit from computation, SPICE on the other 
hand does not exploit the inactivity of the whole lowex bank 
of shift registers and performs the analysis blindly oyer the 
whole circuit in each time step. 
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5*0000+00 
5,0000+00 
5,0000+00 
5*0000+00 
5,0000+00 
f, 0000+00 
1,0000+00 
|,0fOO+fO 
1,0000+00 
5,0000+00 
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Performance of a hand latch 

VOLTAGE WAVEFORMS AT SPECIFIED NODESI 

« f ^ f * f f « :|c « f f ^ f f ^ 3|t f t f f !(( « f f y f 4t 4; « « t t * « 


• I ' ■ ' ■ - . ■ " " ; ■ ■ :■■■■ ^ ' ■ ' ";■ 

fl ' no 


••■■■'•*-M 0 SIMR' 


TiME(N-SEC) 


VOLTAGES (VOLTS 



0?000E+00 


0.100E40I 


0,300E+01 


0.200E+01 


O, 40 OB+Ol 


O,50OE+Ot 


S.OOOE+OO 

olaooiJoi 

0.300E+01 

0.400E4-0t 

o^eooi+oi 

0.800E+01 

0,?0QE+01 

0,100E+02 

O.llOE-fOa 

0,140E+02 
0,150|»02 

o*floi*oi 

0,l80E+02 

psksi 

$ *' 2 1 0 E'fr 02 


0.240E+02 

0.250E#02 

0*2SO|#02 

M28i:8i” 

0*290E#§2 

O,300Eft2 

kinittr 

?:H8i:8i 

hiniin 

0 ^ 390 Bt 02 — 

8:1881:81 

o:j|o|t |5 

0,490l'»-02 

0*SOOE+02 
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A MAND LAfCH 

Table ii-total power versus time 

waveform of total power DISSIPAfEO IN CRT. BLOCK I 

TIMe(N»SEC) total POWER TOTAL POWER CU-WATTS )*•••»-> 


f If 


KQSIKR 




0 , 100 E +03 


0 . 200 E +03 


0 , 3001+03 


0 , 4001+03 


2*f9SE*5l 0,264E+03 

0.600E+01 0,402E+03 

i:§oolJ8| 8:|ofi+?l 

ojup+oa o527ii+ol 

:H itgra ipiisi 

• r 

oh7p+o| l^pel+ol 

0,1801+02 0^276E+03 

kmtU g; ?p:?| 

Ol2|oi+02 0^2761+03 
0,2301+02 04276E+03 

2*1121*92 2J22§®*03 

0,2501+02 0r276E+03 

0,260E+02--0+|76E+03 
0#270E+02 0i|9|E+03 

0,pOE+oa— 0j355C+03 

: ? !:?! 0mil 

: J iJSi 

2«ll2i*22 2-;2§k*2l 
2*1121*22 2»i2lS*2l 

o!l36oi+o2 — oifci?fi+o3- 

g:l?g|Igi Bmxu 

2*^21*22 25lffl*2l 

0,4201+02 0,2761+03 

0,4301+02 0l276i+0l 

2-il2i*2i 2?i2§i*2l 

0 ]* 450 E +02 Qi!i, 276 C+ 0 | 

0,4601+02 0,,276B+03 

0,4701+02 0^2761+03 

mBi lillHI 




'Gi- 

vi'S' 


CiNCLUOfO 




§*♦♦♦♦♦♦2 3 JVug 8 ♦*»♦♦♦♦♦ SPiCf 2G,6 3/15/83 I 

onand latch I -ftS. 

m*** INPUT LISTING TEMPERATURE = 27.000 DEG C | 


» CKT, DESCRIPTION 

.model driver NMOS VTObIV KP*55.2E-6 

.model load NMOS VTO»-4V KPa6.9t-6 

,SU|CKT NAND 1235 

ML 5 3 3 0 LOAD 

MDl 3 1 4 0 DRIVER 

M02 4 20 0 DRIVER 

Cl 1 0 .02PF 

C2 2 0 .02PF 

C3 3 0 ,02PF 

C4 4 0 .02PP 


.ends, 

H 1 4 3 
X2 3 2 4 
♦ inputs 


NAND 

NAND 




VDD,5,,0 DC 5 

Vi i 2 PULSECS ;'5 INS 5NS 5NS 20N5 lOONS) 

V2 2 0 DC 5 

* cdntrol 

.nodeset V(3)a.s V(4)a5 
.tran^ins sons 

* output REQUESTS 

.print TRAN V(i) V(3) V(4) 

.PLOT TRAN V(l) VCh V(4) 

.End I 

Auo 8 ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»#***# SPICE 26*6 

®nand latch 

04»4» MOSPET MODEL PARAMETERS 




TEMPERATURE m. 27,000 DEi C 


DRIVER LOAD 

NMOS NMOS 

1.000 1,000 
1.000 ^- 4.000 

t 2D*05 6»90D«'06 

♦♦♦23 Auo 8 SPICE 26,6 


3/15/83 


TEMPERATURE 


INITIAL TRANSIENT SOLUTION 


RfQDi 


NODE 


NODE V0LTA6E 


NODE VOLTAGE 


node VOLTAGE 


NODE VOLTAGE 




0.0000 


( 1) 5,0000 ( 2) 5.0000 C 3) 0,5357 ( 4) 5.0000 ( 5) 




voltage source currents 
name current 




vdd -s^saoD-os : 

?'H - VI 0 'JoooDtoo ' 'v'; ' ' • 

-tV;.,A''V2 Oi*o6oD+00 "'-''"V/' ' 

optsMiita imfeifii.tMM 








JOB CONCLnOCD 

f: TOTAL JQ^ TIME 6,45 
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table I » ^mLTAGES VERSUS TlWg » 


—-■ — -MOSIMR- 




TIMECN-SEO 


NODE VOLTAGES C VOLTS) 


MODES! 
O.OOOE+OO 
a,30t)E+oo 
a.600E+00 
U.900E+00 
0,1?0E+01 
0,l50Et01 
0.180F+01 
0,210E+01 
0.240E+0i 
0.270E+DI 
O.SOOE+Ol 
O.iBOEROl 
U,360E+01 
0.390E+01 
0,420E+Oi 
0.450E+0t 
0,480E+0l 
0.51dE+0i 
0,S40E+01 
0.57aE+01 
0*600E+01 
0,g30E+Oi 
0,660E+0i 
O.69OE+01 
0»720E+01 
0»75QE+0l 
0.780E+01 
0,81QE+0t 
O.84OE+01 
0,870E+01 
0,900E+0i 
0*930E+0|" 
0,960E+0i 
0*990E+01 
0.102E+D2 
0,105E+02 
0.108E+02 
0*111^+02 
0.114E+02 
0,ll7E+02 
y.l20E+02 
0.123E+02 
0*126E+02 
0,129E+02 
O,i32E+02 
O*'135E+02: ■ 
0.138E+02 
0a41E+O2 
0,144E+02 
O*i47£'f02 
O.fsoE+ol 


J0» CQMCLUDED 


1 

5, 00 (! 

5.000 

5 . 0 0 

5.000 

5.000 

5.0 00 

5.000 
4,91 0 

4 . 0 4 0 
4.370 
4,1 00 
3,830 
3,560 
3.290 

з, 020 
2,750 
2,480 
2,210 
1,940 
1,670 
1,400 
1,130 
0 , B60 
U.590 

и. soo 

0,500 

0,500 

0, 500 

0,500. 

0,500 

0 ,500 

0,500 

O.SOy 

0,500 

0 • 500 . 

0,500 

0,500. 

0,500 

0,500 

0.500 

0,500 

0 , 500 

O.SOO. 

0,500 

0,500 

0,500 

0,500 

0,500 

0,500 

0,500 


2 

0.5 36 
0.536 
0,536 
0.536 
0.536 
0.536 
0.536 
0,536 
0,536 
0,536 
0,5 36 
0,536 
0,536 
■0,.536 
0,536 
0.536 
0,537 
0,541 
0,553 
0.577 

0.714 
0.876 
1,202 
1,764 
2.379 
2.935 
3,405 
3,788 
4,088 
4,276 
4,278 
4,042 
3.566 
2.872 
2.048 
1,401 
1.026 
0,861 
0,835 
0.92 3. 
! ,162 
1,653 
2,258 
2,829 

3,717 

4,033 

4.245 

4,281 

4.082 
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MOS RIMG OSCILLATOR 

voltagf wavefophs at specified nodes; 


IlMECN-SEC) 


voltages c VOLTS 


O.OOOFfOO 


0.1,00E+01 


0,200E+Oi 


^ 'W iPhflj#' 'f 

O.lOOE+00 
0.600E+00 
0,900E+G0 
0,120E+01 
0,t50E+01 
0,180e+Ol 
0.210E+0! 
0,240E+01 
0,270E + 01. 

0 • 300£+0i 
0,330£+01 
0,360E+01 
0,390E+01 
0*420E+01 
0,450Et01 
0*480E+yi 
0*510E+01 
0,540E+01— 
0.570E+01 
O.fOOg+Ol 
0,630E+01 
0*§§0E+0t 
0,690E4-o1 
0,720£+01 
0,7SGEh-01 
0.780E+01 
f,il0S+01 

S:i!Si:g} 

o,9r“ 


0,9|0Et0„ 
O'* 96 get tt 
O.fSOETOl 
0,102E+02 
0.105E+02 
0,108£4.02 
0,lUfe.+02 
0*ll4£+02- 
G*ll7fc-|.02 
0,120E4-02 
0.123E+02 
0.126E+02 
0,129E+02 
0,132E+02 
0.135E+0i 
0,138E+02 
0*14U+02 
0.144e+02 
0,147 Et 02 
0,*l50g+02 


aOi CQlfCtiUOlD 
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table H-THTAL Pn«(ER VEkSUB time 

waveform of total rower niSSTPATFD IM CKT. BLOCK S 

TIMeCN-SEC) total PUWEIR TOTAL POWER (IJ-WATXS) —•—> 


lie 


0 , 000 E+ 0 u 





U,llt»F4D4 


0 • 300E+od 

M, ll0F4R'-i 


0,600i+U0 

0,ltUF404 

# 

0,90RE+00 

y, iiuf:404 

• 

0.120E401 

y.llOFfOA 

fH 

0, 150E+UI 

U, 11 UF4f'4 

m 

0,180E+0I 

U.110F404 

It 

D,?ioi:+oi 

a.ilOF+04 

ft 

0.240E+0t 

U, 11 UF 4 O 4 

# 

0.270S+01 

0,11UF4C4 


0,300l+01 

O.IIOE 4 O 4 

M 

0.330L+01 

C,109F404 

# 

0,360£403 

u,109E404 

ft 

0,390E4Ut 

L,108F404 

ft 

0,4200401 

g,10§F4O4 

# 

0,450E401 

S,107F4G4 

m 

0,480E401 

0 , 1C6F4C/4 

ft 

0.510E401 

0,105F404 

V 

0,540E401 

0,108F404 

ft 

0,570E401 

0,114F404 

ft 

0,600E401 

0,122F404 

' ff' 

0.630E401 

0.132F404 

• 

0,660E401 ■0.139F404 

ft 

0,690E401 

&,137E404 

V 

0*720E401 

0,1208404 

p 

0, 7501401 

0, 1118404 

p 

0i*78OE4Ol 

0,U7E404 

• 

0,8l0g401 

0,1298404 

m 

0.840E401 

O.139E404 

' # 

0,870E40t — 0.i42F404 

m 

0,900E401 

0,129E404 

p 

0,930E401 

0,1218404 

p 

0.960E401 

0,1258404 

p 

0,990E401 

0*135E404 

m 

0, 102E402 

0. 1428404 

m 

0.105E402 

0, 136F4C4 

p 

0,108E402 

0,119F404 

m 

0,U1E402 

0,112E404 

p 

0,ll4a402 

0,1178404 

t 

0.tl7E402 

0,1268404 

' P 

0,1208402 

0,135S404 

' t 

0,12 3£402 

0,1378404 

p 

0, 12 6 E 40 2 

0,1258404 

f 

0,129E402 

e,tl4E404 

P 

0,132E402 

0, 1165404 

P 

0.135E402 

0*1^65404 

f' 

0.138E402 

01 1375404 

9 

0,141E402 

0.1425404 

p 

0,144£402 

0,1315.404 

p 

0,147E402 

0, 1215404 

f. 

0,1506402 

0,124E404 

. P 


JOB CONCLUDED 





0,ia8E+04 


0.360E+03 


0.720E+03 


0,144E+04 


0,180E+04 


To-J-cl^ jsotoe^-’ 

JU ssi p>aJt^a^ ’ 
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KIMns KING DSf'IT.r.^ATnH 
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TTMP 

X 

VCD 

0(7) 

O.-UOE + OO 


A,35Pf>0l 

3,''00e-1^‘ 


5 • 35PE— 1‘ 1 

6,':0CF>1P 

5. '■if O') 

35RE«(^1 

q.^^ooe-io 


A,350e-61 

l,?0PE-0« 

5 , r 

5 , ISPE- jt 

1 .50PE-0'’ 

, 0(}f">t:;+ l""i 

5, 35PE-OI 

t .fO^E-OO 

5,f''0OE + 0 0 

5.35OE-DI 

2 ,ioof>o 2 

I.P'lf’L+UO 

5,35PE-0t 

240OE-O2 

1 . A40L4-0n 

*'.3S9e-0i 

^.TqOE-OP 

4 , '»70i''>00 

5, 359e-0l 

3,f;00E-0P 

4,1 O^-’K+uo 

5,159E-01 

3.300E-OP 

3,R3fiE+0'’ 

5, 35RE-OI 

3,6oOE-OR 

3,560^4-00 

5,35«E-01. 

3,900E-00 

^.?90E+U0 

5.350E-01 

4 , 20DE"0® 

1,0201::4.00 

5,150E-01 

*4 , 500E*09 

2.*'50t:+oo 

5.367E-01 

^.foOE-O^ 

7.48r)E+C'' 

S,38(>E-01 

5.t0DE-0R 

7,7loE+00 

A,441e-oi 

5,4oOE-OP 

1 ,Q4'*y+0O 

5,6o5E-01 

5,70Oe-09 

1 .670E+00 

5,07tE-0l 

e.^’OOE-OR 

1 , 40OE+00 

A.716E-01 

6,300E-09 

1 . 1 30E4-0O 

7, 7 J RE-01 

6,600E-09 

0, AO OS-01 

t ,150E+00 

6.900E-OR 

A,00OS-0i 

t ,73tS + 00 

7.200E-09 

5,00''E-01 

7,4512+00 

■ 2.500E-0R 

5.oooe-oi 

3.098E+U0 

T.fOOE-O*^ 

<1,oonL-01 

3,630S+O0 

8,100E-02 

5.O0OE-0I 

4,032E+0O 

?*.400E*09 

5,00OE-0t 

4,274E+oO 

P*I00E«<}9 

A, oOOE-oi 

4,2oRe+00 

R.COOE-OR 

5,00r'E-01 

3.932E+0O 

9*300E*^09 

5,OOOE-0! 

3,749e+00 

P.600E-yP 

5.O0OE-01 

7,335E+0O 

9,900E-09 

5,oook-oi 

1 ,493E+00 

i ,©20E-0R 

5,oooe-oi 

R,R33e-0l 

l.fS0E-'0« 

5,oont>oi 

7i,‘360E-Ol 

l,fl0E-OR 

5.000E-0! 

7, t40E-0l 


5,O00E-01 

1 ,007E+00 


5,OO0E-O1 

i,445E+oO 

; l*t70E-08 

5,000E-01 

7,09RE+0O 

! , 200E*0R 

5.000E-01 

7.787E+00 

1,230E-0R 

S,000G-01 

3,380E+00 

1 ,260E-09 

5,00OE~O’ 

3,R4RE+00 

t .290E-OR 

5,oooe-c’i 

4,177E.*-00 

1,320E-0« 

5,O0OE*(M 

4,3o4E+00 

l,350E-0fi 

5,0f>0E-U1 

4,1 33E+00 

1 ,380e--OR 

5,0UOE-01. 

3,6l5E+00 

1,410E-0*=* 

R,ooos«(^!^ 

7,78lE+00 

1>40E-0P 

5,(in»oe-r-i. 

1 ,85fiE+0O 

l,470E-0fi 

5,O0OS-0t 

l,176E+00 

1,500E-08 

5,O(»OE-0l 

R,253E-01 
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^ILFGFND: 

*t V(1) 

+! VC ?3 
X 

ttmf 

o.<‘^on+Du 

3.i.n0D-10 

6. < OD-IO 
9.^>*'OD-10 
i,2’^02«09 
1.5202*^9 
l,8*'0D-09 
|.lf'0D-O9 
2.4«^on-09 
2.7^0D-09 

3,5;^'0D-D9 

3,3f‘0D-.09 
S.6'*0D-09 
3.9'’*0D-09 
4.2^on-09 
4.5^00-09 
4,al»0D-09 
5.1^>OD-09 
5,4non»09 
5,7‘70D-09 
6*eF0D-09 
6, 3000-09 
6,6£0D-09 
6,9»*00-09 
7.2^00-09 

7. §£00-09 
7,8<SOO-09 
8.1000-09 
8.4000-09 

8 . feS - 2 §- 


9.#»0D-^O9 

9.||gg-09 

JM'oo-os 

“^•11)500-08 
1 .©800—08 
1.1100-08 
1,1400-08 
1,1700-08 
l,2O00-08 
1.2300-08 
1.2600-08 
1.2900-08 
1,3200-08- 
1,3500-08 
i.35{)D-08 
i, 4100-08 
1.4400-08 
1 . 4700— 0 8 
l,5f00-08 


vm 


0,0000+00 


2,OOOD+O0 


, 0000+00 


6,0000+00 


8,0000+00 


5,0C-fM) + 0f' 
'5,00‘''D + IT 
5.00^0+00 
5,ry(;U + 0O 
B,oeou+uo 
5,OOOD+UD 
5,OOPD+Or 

4. R1CD+(!0 

4.6400 + 0'^ 
4,370D+0D 
4.10OD+t)C 
3.8 3fiO*OP 
1.56OD+0f' 
B.^gf^o+up 
1,020D+0C 
7,750D+0P 
2.^800 + u^* , 

7, '>100 + 00 
1,0400+00 
1,6700+00 , 

t,4U0D+UO 
1,1300+00 
8,6000-01 , 
5,0000-01 
5,0000-01 
5,0000-01 
5,0000-01 

5, nu-^o-oi 

5,0000-01 

■5,0000-01 

5, 0000-01 

5.0CO0-01 

5,0000-01 

5,0000-01 

5,0000-01 

5 , 0000-01 

5,0000*01 

5,0000-01 

5,0000-01 

5,0000-01 

5,''000-Ol 

5,0U00-0l , 

5,0000-01 

5,O00i)-Ul , 

"5,0 000-01 — 

5,Oim!)-01 
5,0u0u-01 , 

5,0000-01 , 

5,OUOi)-01 
5,0000-01 , 

5,0000-01 . 
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JOB CONCKOOEO 
TOTAL 108 TIME 


6,06 


2 . 4 . 3 e . 






f 




120 
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VOLTACIE waveforms at SPFCTFTEn NOHEf?; 




TTMECM-SED 


vot.tages f vnijTs ) — — 


MnSTMR 


? 10 I J V nr V A,l 


V * 4 - V r T 4 


> m r I 


M • W X 4. 



0 ,000E+00 

o»iooF+oi 

0,200F+(31 

0,300F+01 

0,400E+01 

0,500E+01r— 

0.600E+01 

0,700E+0i 

0,800E+01 

o.gooc’+oi 
0,l00E+02 
0,li0R+02 
0.1201+02 
0*1 101+02 
0.1401+02 

0,1701+02 
0,1801+02 
0, 1801+02 
0,2001+02 
0,2101+02 
0,2201+02 
0.2701+02 
0,2401+02 
0*2501+02 
f*2§0E+02 

in |7 01+ of — 
|*|f 01+02 — 
|*|OOE+O| 

0*3201+02 

0,3301+02— 

0,3401+02 

0.3501+02 

0,3501+02 

0.3701+02 

0,380P’+02 

0 , 3501+02 

0 , 4001+02 

0.4101+02 

0.4201+02 

0,4301+02 

0,4401+02 

0*4501+02 

0,4601+02 

0*4701+02 

0,4801+02 

0 * 4501+02 

0.5001+02 












CONCliUDFD 





vet i 

V()? J 

veil) 

VC9) 


TTHf? 


^ "f* 5* t * 


tff n 




$<»)' 

0 . y«Uft+00 

1., /’on«09 
i!.4'*0n-09 

3. at^on-09 

4. -inoR-09 

5. ^:i0n-09 

6. (^non _ 

7.. <.^^0D-O9 
B,',/?^aO«09 
9.vft0D-09 

1 , M«on-08 

l.i«'»on-o8 

1. 2^100-08 
350n-O8 
4«OD-08 
S%^0D-08 

(j 01)0 ‘•08 

1,7000-08 
i.8«0n-08 
I , 9O0D-08 

2. «s?0!)«0t 
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CHAPTER 5 


CONCLUSION 

5 . 1 ADVANTAGES ACCRUED BY IviOSM: 

MOSIMR as a special purpose program to perform time 
domain analysis of NMOS digital circuits has been able to 
gain several advantages., 

5.1.1 Improved Speed of Analysis: 

As apparent from the results in Chapter 4 MOSlMR 
proves its efficency in the analysis of comparitiveiy large 
circuits. Node partitioning and effective exploitation of 
circuit latency are chiefly responsible for it being more 
than ten times faster than SPICE 2G even in the case of medium 
sized circuits (such as the two-bit Johnson ring counter). 
Single iteration of GS and NR employed by MOSIMR makes the 
program faster than other special purpose iterative simulators 

5.1.2 Data Storage Requirement, Minimal: 

The program MOSIMR uses equational device models, thus 
avoiding the enormous storage required for tabular models [lO] 

Using nonlinear Gauss-Seidel method the program only 
needs to maintain a one-dimensional voltage array independent 
of simulation period. Vi/hile RELAX [ll], [l2j, using waveform 
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relaxation method, needs to store complete waveforms of N 
circuit nodes over a whole simulation period (consisting of 
T time steps) in a array. Thus memory requirement in 
vyaveform relaxation method is proportional to the period of 
simulation. Estimated storage needed in MOSIMR fpr a large 
circuit with 1000 nodes, is 55K words of core, independent 
of simulation period. This may be compared with the data 
space requirement of 2,4M bytes by RELAX to analyse an MCS 
circuit with 1000 nodes, over ICX) time steps [llj, 

5.1,3 Multifarious Applications of MOSli®; 

Though a special purpose simulator, different features 
of MOSIMR can be utilised to study various circuit aspects as 
followst 

i) MOSIMR computes the initial condition of the circuit 
under steady state, based on a full-fledged logic simulation. 
Hence it can potentially be used to simulate logic of a 
digital NMOS circuit under any steady state condition. Logic 
function of a complex combinatorial circuit can thus be found 
out with the help of MOSIMR applying different input combi- 
nations, Along side, quescent power dissipation of a circuit 
under any steady state condition can also be obtained, 

ii) Speed and power dissipation becomes important factors 
in the design of sequential circuits and finite state machine 
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structures. The usual timing analysis feature of MOSIMR 
can be used to study response of such circuits with clocking. 
This way the maximum allovaable clock frequency can be found 
out for a design, 

iii) The external definition of input (EXDEF) allows 
arbitrary waveforms to be applied to the circuit. Taking 
advantage of this feature a large circuit can be brgken up 
into blocks which can be analysed independently, over the 
entire simulation period. The output from the previous 
block is printed in the output file OUT.I'CCT. This output goes 
as EXDEF input to the next block and OUT .NXT is renamed 
INP.PRE to act as an input file during the analysis of next 
block. Data, read-in from INP.PRE, in every simulation step, 
defines the EXDEF input of the next block, 

iv) MOSIMR can detect the signal-wave profile, as it 
travels through a large circuit. Hence a novel application 
of MOSII\® would be as a signal tracer in a large circuit. 

v) Interactive user control over ihe simulator allows 
the user at the end of a specified simulation period to 
decide whether to STOP, RESTART or COi'HIMJE furthej? simulation. 
This feature can be used in dividing a long simulation period 
into subintervals. As simulation is completed smoothly in one 
subinterval, the user can direct the simulator to continue 
to the next. 
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5.2 LIMITATIONS IN MOSLMR AIO SUGCESTIONS FOR BETTERi.ENT: 

Circuit initialisation related with the *On* pass 
transistor makes an approximation in computing the source 
output voltage, when the drain input voltage is ujiknown and 
the input is at logic *0* (Sec. 3,3,2). This may be avoided 
by conducting the voltage initialisation trace alortg the 
signal flow path, as done in logic simulation. 

Handling of floating capacitance in MOSIMR puts a 
constraint over the time step h (Sec. 2.2.]j). Better way to 
handle these v/ould be as done in RELAX [ll]. 

In MOSIMR the link graph structure of the circuit is 
built up with the help of linked list. Implementation of the 
linked list data structure would have been easier with PAiXIIAL. 
Also recursion in PASCAL v/ould support multiple-level 
subcircuit nesting. 

Alternative to the linked list representation of the 
directed circuit graph, is the dependency matrix approach 
(Sec. A.'i-n.S) . In this method the program can be easily 
adapted to simulate circuit in transistor level and allow 
definition and use of complex logic gates. 

In MOSIMR accuracy and stability are jeoparAiz-ed with 
large value of simulation time step h. If stability is our 

\ 
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concern a method should be formulated to redefine GG (Sec, 2, 4 ) 
in a way such that single iteration of Newton Raphson (NR) 
gives an asymptotically correct result, with h— On the 
other hand to ensure accuracy i€i iterations must be carried 
to convergence in each Gauss Seidel (GS) relaxation. Also 
several iterations of relaxation must be performed per 
time step, to gain at least partial convergence of solution, 
in the analysis of non-unidirectional circuits, 

5,3 SCOPE FOR FUTURE VVORKs 

MOSl]\iR has been built to simulate NMOS circuits with 
basic ratioed logic gates. Further effort can be directed in 
expanding the simulator for ratioless gates, complex gates 
and finally for CMOS circuits. 

This will involve mainly incorporation of transistor 
level simulation (over and above basic gate-level simulation), 
introduction of new device models and definition of new 


circuit elements. 



APPENDIX A 


ELUCIDATIOSI of few TEFU^iS AND CONCEPTS 

Here, some terms and concepts, often used and referred 
in the preceding chapters, have been elucidated with defi- 
nition, brief explanation/ derivation and illustrative 
e xample s • 

A,1 OVERHEADS: 

By 'overheads* is meant here the extra processing (not 
directly connected v\?ith bare circuit analysis), which is 
required in order to make a program general purpose arjd user 
oriented. 

It may imply extra CPU time as well as extra storage 
for simulation. 

i) A language translation stage at the beginning of a 
simulator program to allow the user to specify input 
data in free format j 

ii) Storage of models of a large variety of devices to make 
the simulator generally useful; 

iii) Extra accuracy in device models, than is often required 

iv) Elaborate outputting - giving details in listing, in 
which the user is often not directly interested; 



135 


t are, to name a few, which principally cause overheads in a 
general purpose circuit simulator, 

A. 2 MCROMODEL: 

A circuit macromodel is a unit within a circuit which 
is actually a physical conglomerate of several el^nt^nts and 
devices. So it is rather a complex element which dan be 
modeled either on the basis of component element models (i.e, 
device level macro-modeling ) or simply by characterising the 
input/output functional relationship without regard to the 
physical interior of the unit (i.e, functional macromodeiing) • 

In mOSIMR a macromodel is termed, a circuit element, 

A, 3 RELAXATION METHOD: 

General scheme of solving a system of algebraic 
differential equations by relaxation method [ll]f El^] 
consists of two major processes, namely the assignment parti- 
tioning process and the relaxation-iteration process. 

F(y,y,u) =0 (A, la) 

E(y(o)-y^) = O (A. lb) 

(A. la) is the system of equations to be solved with (A. lb) 
as initial conditions, where 

y(t)Cf^ is the vector of p unknown variables at 
time t. 
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is the time derivative of y at time t 
and u(t)CR^ is the vector of r inputs at time t 

A. 3,1 Assignment Partitioning Process: 

In assignment partitioning process, each unknown 
variable is assigned to an equation of (A, la) in which it 
is involved., Thus (A. la) is partitioned into p disjoint 
subsystems: 

Fl(yi» yj^» d^, u) = 0 

'^ 2 * ^2’ u) = 0 (A.2a) 


^p’ '^p' “) = 0 

E(y(o) - y^) = 0 (A.2b) 

where yjj_(t) is the variable to be solved from 

= 0 (A.2c) 

with d.O<V^, y 2 ,...y^_p 

Consisting of other variables present in eqn, (A, 2c) 
which are considered as inputs to the subsystem, represented 
by eqn. (A. 2c), with their trial values. 

Consideration of d^ helps to decouple the eqn. (A. la) 
into those in eqn. (A. 2a) and hence the variables included 
in d^ are called the decoupling inputs of eqn. (A*2c). 
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A, 3, 2 Relaxation-iteration Process: 

The two most comaonly used types of relaxation methods 
namely the Gauss Seidel (GS) and the Gauss Jacobi (GJ) are to 
be considered here* 

The initial trial solution for y(t) is a guess* During 
each iteration each decomposed subsystem of eqn. {A*2a) is 
solved for its assigned variable using the trial values for 
the variable and the corresponding decoupling inputs* 

For the GS relaxation, the solution obtained for a 
variable by solving a decomposed subsystem is immediately used 
to update the trial values - so that it can be used in the 
solution of other decomposed subsystems which follow in the 
same iteration step. 

For the GJ relaxation, all trial values are updated 
only at the beginning of next iteration. 

The relaxation process is carried out repeatedly as 
the trial values converge to the actual solution. 

The convergence criteria for relaxation-iteration are 
given in [llj and these require diagonal dominance of the 
Jacobian of the discretised equations obtained frofti 
eqn , ( A.la) . 



138 


A *3. 3 GS/GJ Time Advancement Formula; 

This is a recursive formula for solving (A. la) in 
successive time steps Vi^ith only one relaxation iteration 
performed per time step. As this is a kind of numerical 
integration the method is called Guass-Seidel/Gauss-Jacobi 
Integration algorithms, according as the relaxation scheme 
used is GS or GJ, 

Illustrating the time advancement scheme for a linear 
system described by linear state eqns,; 

X = Ax (A. 3a) 

and x(o) = x (A, 3b) 

where x(t)CR'^ the state vector to be solved at time 

t consisting of n variables, 
the coefficient matrix = L+D+U 

where L is strictly lower triangular 
D is diagonal 

and U is strictly upper triangular 
The G.S. integration formula is - 

[ I - h(D+L)3 = [l + hU] x^ (A, 4a) 

and the G-J integration formula is - 


ll-m] x^^3_ = [l+h(L+U)3 


(A. 4b) 
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where I is the identity matrix, 

^t+1 state solution at t+1, 

is the trial solution, which rather being solution of the 
previous iteration is the solution of previous time step (as 
only one iteration of G-S/G~J relaxation is performed per 
time step), 

A, 3, 4 Network Ordering and its Effect on Timing Errors: 

In G-S/G-J integration algorithm (Sec. A. 3.3} only 
one iteration of relaxation is performed. Now we should 
recognise a fact here and that is circuits with unidirectional 
properties can be analysed exactly by performing only one 
iteration of G-S relaxation if the computation of the node 
voltages are ordered in such a way that the decoupling inputs 
to each equation of the partitioned system axe evaluated, for 
the present tine step before the particular equation is solved. 

This ordering of computation in the case of an unidire- 
ctional circuit, for which a link graph can be built (A. 4. 4) 
(where computation of a node depends on the higher order nodes 
alone and not on the lower order nodes (A«3,9.3)) is naturally 
according to the precedence order of signal flow. 

Thus network ordering has a profound effect on the rate 
of convergence of G-S method. However, or dering of computation 
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does not have any effect on G-J method where the decoupling 
inputs retain the satoe values throughout the whole iteration 
step. 

Now timing error is a phenomenon that is always 
evident in G-J method and becomes evident in G-S case vdien 
computations are done in an arbitrary order or even with 
network ordering when the circuit is not unidirectional. 



Fig. A.lb: Circuit with 
feedback 


Fig, A.ld: Circuit with floa- 
ting capacitor 


Considering the example of unidirectional circuit in 
Fig, A, la the ordered G-S would proceed to compute the nodes 
in the sequence 1-2— 3-4 following the flow of signal in the 


circuit 
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at 6 \ ^ timing error of one time step is 

incurred. 

Similarly computation of 2 at depends on 3 

tj^ in Fig, A.lc and Fig, A, Id, incurring timing errors 
of 1 simulation step in each case, 

A*3v5 Latency: 

For most circuits the fraction of nodes which change 
their voltage values at a given point in time , decreases as 
the circuit size increases. For circuits containing over 
500 MOSFET^s, fewer than 20A [l2] of the node voltages 
change (significantly) over a simulation step. This ina- 
ctivity in circuits have been termed ’latency'. 

When circuit analysis is scheduled in such a way that 
at any time only the active nodes are evaluated, leaving out 
the redundant calculations of the inactive parts, we tall 
it 'exploitation of latency', 

A. 3. 6 Leading Element: 

In the linked structure based on the concept of 
observability (Sec. A, 3, 9), node N in Fig. A,2 is observed 
by all the three elements 1,2,3 which together form a chain 
where 2 may be thought of as the next brother (N3) element 
of i;3, the NB of 2 and so on. The first element 1 v^hich 
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is recognised as the successor (succ.) element at M,is 
termed the leading element of the MB chain from which a 
trace scheduler can visit the subsequent 1'JB‘s following the 
NB pointers. 

Similarly in the fan-out table representation of 
dependency, the first element, in the list of fan-out 
elements, recorded for node N, is the leading elaioent. 



» 

t 

i 


Leading element among the 
fan-out elements at N,is the 
one topping the fan-out list, 
namely. the pass transistor 
marked *1* here. 

Fig. A. 2. a Fig, A.2.b 

A*3,7 Selective Trace: 

Whenever the voltage at a node changes, it is possible 
to schedule all the elements observing the node to be processed. 


Leading element in the chain 
of next brother elements, 
following gate G at node N, 
is the one marked •!' 
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In this way the effect of a change at the input of a circuit 
may be traced as it propagates to other circuit nodes via a 
linked list/fan-out table etc* 

Since the only nodes processed in this way are those 
which are affected directly by the change, this circuit trace 
procedure is selective and is known as ’selective trace*. 

A, 3, 8 Scheduler: 

A scheduler is a program which plans which nodes are 
to be computed and in what sequence in a particular time 
step of simulation. 

A. 3. 9 Controllability and Observability of Nodes: 



A piece of network showing the 
node of interest N and its 
neighbouring nodes N1-N5 which 
have controllability and/or 
observability relations with N. 



denotes control re- 
lation. 



denotes observabi- 
lity relation* 


Directed graph represen- 
tation of the piece of 


network. 
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Voltage at the test node N in Fig. A, 3a can be changed by 
varying the voltage(s) at any. ox more of the nodes ni, N2, 

N3 and N4, whereas N5 does not have any effect on N« So we 
say Nl, N2, N3 and ^4 are the nodes vdiich contr ol N. 

On the other hand any change of voltage at node N 
cannot be felt at Nl, N2 and N4 but it can surely be felt 
at N3 and N5« So N3 and N5 are said to be the nodes which 
observe N* 

These relations are conveniently represented by a 
directed graph where the nodes of the graph correspond to the 
circuit nodes and the edges of the graph correspond to the 
controllability and/or observability relations between nodes 
through elements. The directed graph-equivalent of the piece 
of network in Fig* A. 3a is given in Fig* A. 3b, 

A*3.9*l Reachability ; 

A node j is reachable (by controllability/observability 
relation) from node i if and only if there is at least one 
directed edge from i to J [l8]* 

A. 3, 9*2 Strongly connected graph : 

A graph is strongly connected if and only if every 
node is reachable (Sec. A.3,9.1) from any other node. This 
means that in a strongly connected graph, every node is in 
at least one loop [l8]* 
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A*3,9i3 Maximal strongly connected subgraph : 

A maximal strongly connected (MSC) subgraph is a set 
of all possible nodes that are strongly connected (Sec, 

A. 3, 9*2) with each other, in a graph. 

This designates the maximum set of nodes that are- 
involved in conflict (Sec. A*3,9*4) from controllability and 
observability view point [9]* 

A,3,9*4 Order of nodes 

Higher order nodes are those evaluated earlier than 
the lower order nodes [9], 

A, 3, 9,5 Conflict : 

Conflict is said to .exist when computation of a higher 
order node depends on any lower order node(s) [9], 

A, 3, 9 ,6 Link graph ; 

A link graph is one that contains no strongly connected 
(Sec, A. 3, 9. 2) or unconnected subgraphs in it. It is there- 
fore an acyclic digraph, (i,e, loop free) and the nodes in it 
can be arranged in an ordered list. Physically this corres- 
ponds to a unidirectional circuit where in principle, with 
ordering of elements, nodes can be computed without any 
conflict (Sec, A. 3, 9, 5). 
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A ,3, 10 Resolution of MSG*S; 

MSC's (Sec, A»3.9*3) which cause conflict due to 
presence of loops, must be resolved by breaking the loops 
(known as decyclization in graph theory) arwi introducing 
additional control and/or monitor points associated with each 
broken edge - to maintain equivalent of the derived graph. 

In fact all edges entering the node of conflict are 
broken this way and additional control and/or monitor points 7 
introduced. 

The node of conflict is hereby partitioned and a 
subsystem is formed which includes the node in question and 
ail the Control points introduced. In such a subsystem node 
N can be evaluated independently with the help of higher order 
nodes and the additional control points introduced in the 
subsystem block. 

Hence node partitioning can resolve MSC's and allow 
a link graph structure (i.e, an acyclic connection of the 
above subsystems) to be formed, 

A, 3, 11 Node Partitioning and Ordering: 

Node partitioning and ordering can be dong wholly on 
the basis, of controllability and observability relations 
respectively, of nodes. 
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A subsystem obtained after partition should contain 
the node of interest and all the controls on that node i.e. 
all the fan-in elements (Sec. A. 3.11.1), 

Ordering of partitioned subsystems determines the 
order in which node(s) to be visited after computing a node i. 
Naturally those nodes which observe node i are the ones to 
be visited after updating node i, 

A linked list [l9] impl^entation of the ordering is 
possible where the pointers depict observation edges that lead 
to the next node(s) to be visited, 

A fan-out list for each node is an alternative option 
(Seq, A.3.11..2). 

In Fig, A. 3a for example the node of conflict M, when 
partitioned., gives a subsystem that includes the node in 
question i.e, N, the fan-in elements at N i.e. the NAND gate 
(NG) and the pass transistor (PT) and fan-out elements at N 
i.e, the pass transistor (PT) and the inverter (IG). The 
fan-in, fan-out information for each node can be kept in 
tables or other wise can be elegantly represented with the 
help of dependency matrix (Sec, A. 3. 11, 3), 

A. 3. 11.1 Fan-in Element: 

Fan— in elements of node i are defined as those which 
play Some part in determining the voltage at node i (i.e. 
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those elements which cause some entries of column i of 
dependency matrix to be l), 

A. 3 •11, 2 Fan-out Element ; 

Fan-out elements of node i are the ones whose operating 
Conditions are directly influenced by the voltage at node i 
(i«e, those elements which cause some entries of row i of 
dependency matrix to be l), 

A.3*ll,3 dependency matrix ; 

A signal flow graph which determines the network 
ordering for the processing of nodes, can be generated by 
means of dependency matrix. This is the adjacency matrix, 
of the circuit digraph, of order nxnj where n = total number 
of nodes. 


Entry in the matrix = 1 if there is a directed 


O, When 


edge that joins node x^ to node x^. Otherwise P^^ 

P^j =1 the physical inplication is, node voltage v^ influ- 
ences node voltage Vj via a device equation. That is is i 
fan-in node of x. and conversely x. is a fan-out node of x.. 
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APPENDIX B 

MOSIMH : USER'S GUIDE 

MOSlArlR is a special purpose circuit simulator for tin^ 
domain analysis of digital i'MOS circuits. This manual 
elaborates how to interact with the simulator to get desired 
results. It also points out the constraints and idio- 
syncrasies of the program in its present version. 

B.l. THE INPUT /OUTPUT FILES: 

Input data is read from three files, namely: 

1) CKT.DES - This contains the linked description of the 

circuit to be analysed, in terms of elements 
(table B.l) and subcircuits. 

2) INIT.CQ^ - This contains the initial condition of the 

circuit and the output requests. 

3) INP.PRE - This contains the EXDEF( table B.2) inputs in 

a specific order, tabulated for every simu- 
lation time step, after t=0, over the entire 
simulation period. 

R.e suits are printed in the foliov/ing output files? 

1) VTABL.DAT - This contains the requested node voltages 

printed against time over the whole simulation 
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2) VCURV.DAT 


3) PTABL.DAT 


4) PCURV.DAT 


5) TPCUEV.DAT 


period. Default printing/plotting step is 
around 1 nsec. At most 10 voltages can be 
requested to be printed and plotted. 

This contains the plot of requested node 
voltages against time. Numerals 0-9 are used 
to represent the (l-io) consecutive requested 
node voltages. 

This contains the requested elem. oowers 
printed against time over the entire simulation 
period. Default printing/plotting step is 
around 1 nsec. At most 10 elements pan be 
requested for their powers to be printed and 
plotted. 

This contains the plot of element-power curves. 
Numerals 0-9 are used to represent the powers 
of the consecutive (1-10) requested elements. 

This contains the print and plot of the total 
power dissipation of the circuit, over the 
entire simulation period. The letter 'T* is 
used to plot the power waveform. Default 
printing and plotting step is around Insec. 
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B.2 INPUT FILE FOFilvL^TS: 

The data provided in the input files CKT.DES, INIT.CON 
and INP.PRE, are in alphanumeric, integer or real modes, with 
following formats. 

B.2«l General Formatting of Data in CKT.DES: 

Each line of data constitutes a data card of the input 
deck. The sequence of data cards in CKT.DES is as follows; 

1) Title card - This contains an alphanumeric string of 

maximum 125 characters. 

FORkiAT (3X, 25 A5) 

2) jViAXNOD, NSUB - This contains the total number of global 

nodes in the circuit description (xViAXNCO) 
and the no. of subcircuits (NSUB) descri- 
bed. 

F0RMAT(3X, 2(13, IX) 

3) Subcircuit name - This contains the name of the subcircuit 

to be described. 

FCW/iAT (3X, A5) 

Skip this card and the subsequent subckt. 
description cards (upto FINIS) if i'SUB 
is specified ’0* in the previous card. 
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4) NAi\E, NO - This contains the name (NAi£) and number 

(local) (no) of the subckt. elements 
being specified, 

FORJviAT (3X, A5, IX, 13) 

5) Subckt. element specification ~ This contains the element 

specifications in the following sequence: 

Input node(s) (local), output node(s) 
(local) name and no. (local) of next brother 
(i'B) element(s), name and no, of successor 
(succ.) element(s). 

The exact format in which the above data are 
specified depends on the type of the elero. 
and can be obtained from Table B.l. 

After this elem, specif ication the user 
specifies all the remaining subckt, elements 
in the san^ way, using cards 4 and 5. 

6) FINIS - This contains the reserved word ’FINIS', 

indicating the end of subckt. description. 
FCRfvlAT (3X, A5) 

After this the user describes the remaining 
subckts, (if NSUB > l) in the same way, 
using cards 3-6, 

After all the subckts. have been described 

the main circuit elements are specified 
using cards 7 
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7) NAivlE, NO, LIN(l), 1-1,5 - This card contains the naae 

(NAME) and number (global) (no) a main ckt, 
corresponding 

elem* / to the input and output nodes of 
the element there is a variable (LIN) to be 
specified. Specify LIN as ‘L’ in case the 
element is to be linked with a subckt, elem. 
at the corresponding l/o node. Otherwise 
specify LIN as ’N’, 

FCBMAT(3X, AS, IX, 13, IX, 5(A1,1X)) 

8) Main ckt, elem, specification ~ This contains the element 

specifications in the following sequencej 
Input node(s) (global), output node(s) 
(global), name and number (global) of the 
NB element(s), rame and number (global) of 
the succ, element(s). 

Exact format depends on the elem. type and 
is given in Table B,l. 

After this, the user specifies the 
remaining main circuit elements and inputs, 
in the same way, using cards 7 and 8. 

9) Subcircuit call cards «- Calls to the previously described 

subcircuits are placed after all the main 
circuit elements and inputs are specified. 



155 


Cards 9(a) - 9(d) are used to place repeated 
calls to a particular subckt. 

9a) CALSB, NO - This card contains the reserved wrd *CALSB* 

indicating subckt, call, and the number of 
calls (NO) to the same subckt, FCHMAT 
(3X, A5, IX, 13), Skip the cards 9(a) - 9(d) 
if NSUB in card 2 is 0. 

9b) NAivlSUB,. NOEX - This card contains the name of the subckt. 

(NAMSUB) and the no, of external nodes of 
the subckt, (NC£X), 

FORMAT (3X, A5, IX, 13) 

9c) ((EXTERN(I,J), J=l,3), 1=1, MX) - This card contains 

dummy external nodes (EXTERN(I,1) ) of the 
subckt, and names (EXTERN(I,2)) and local 
nos, (EXTERN( 1,3) ) of the subckt. element 
to be linked with the main ckt. environment, 
at the duramy ext. node, 

FORMAT (3X, 9(13, IX, A5, IX, 13, IX)) 

9d) (EXTNQD (j), J=i,NOEX) - This card contains the set of 

actual nodes that replace the dummy external 
nodes of the subckt. The actual nodes are 
listed in the same order as the dummy nodes 
in card 9(c) 

F0RMAT(3X, 9 (13, IX)) 
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If No. in card 9(a), is greater than 1, then 
the same subckt. is called more than once in 
different main circuit environments. For the rema- 
ining cases the set of actual nodes are listed 
in the same way, after this, using data card 9(d), 

After all ’NO* sets of actual nodes have been 
given like this, calls to the remaining subckts 
(if NSUB in card 2 is > 1) are placed in the same 
way using cards 9(a) - 9(d). 

10 ) CLOSE - This card comes after all the subckt, calls have 
been placed. It contains the reserved word 
'CLOSE', indicating the end of the circuit descrip- 
tion file, CKT.DES. 

FCRMT {3X, A5) 

B.2.2 General Formatting of Data in INIT.CON: 

1) NNODE, WM, WP - This card contains the no. of nodes for 
which path-length stray capacitance would be 
computed (NNODE). 

Alongside, the widths of the metal lines (Mw) and 
poly lines (VVP) are also mentioned, if these are 
different from default values. 

Default value for VVM = 9.0 pm 
and for WP = 6.0 pm 
F0RimT(3X, 13, 2(1X, F3.1)) 
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2) NODE, NTYP, PATHL - This card contains the no. of the 

nodes whose capacitance is to be computed (NOEE), 
type of the node (NTYP) which is either 'm* or 
*P*, indicating metal or polysilicon respectively 
and the path length (PATHL) in pm. (PATHL cannot 
exceed 1000 pm). 

If PATHL is given as zero the default value 
assumed is 30 pm. 

If NTYP is left a blank, default type 
assumed is ’M* . 

Thus default value of stray capacitance of 
any node (unspecified) is ,0081 pF. 

F(BimT(3X, 13, IX, Al, IX, F4.1 
Skip this card and go to card 3 if NNODE in card 
1 is equal to zero. 

If NNCDE > 1 the remaining nodes are 
specified in the same ;vay, after this, using 
card 2, 

3) NO - This card contains the no, of nodes whose 

initial logic states are specified (NO). 

FORiVAT (3X, 13) 

4) NODE, LOGVAL - This card contains the node no. (NODE) and 

its specified logic state (LOGVAL) (either a 
•0' or '1’ )• 
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2) NODE, NTYP, PATHL - This card contains the no. of the 

nodes whose capacitance is to be computed (NOIE), 
type of the node (NTYP) which is either ’m* or 
’P*, indicating metal or polysilicon respectively 
and the path length (PATHL) in pm. (PATHL cannot 
exceed 1000 pm). 

If PATHL is given as zero the default value 
assumed is 30 pm. 

If NTYP is left a blank, default type 
assumed is 'M* . 

Thus default value of stray capacitance of 
any node (unspecified) is .0081 pF, 

F0El£\AT(3X, 13, IX, Al, IX, F4.1 
Skip this card and go to card 3 if NN0E£ in card 
1 is equal to zero. 

If NN(X)E > 1 the remaining nodes are 
specified in the same v/ay, after this, using 
card 2. 

3) NO ~ This card contains the no. of nodes whose 

initial logic states are specified (NO). 

FOPJvlAT (3X, 13) 

4) NODE, LOGVAL - This card contains the node no. (NODE) and 

its specified logic state (LOGVAL) (either a 
'0' or *1’ )* 
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FCEiviAT (3X, 13, IX, II) 

Skip this card if NO in card 3 is equal 
to zero. 

If NO > 1 the remaining nodes are specified 
in the same way, after this, using card 4. 

5) nodint — This card contains the no. of nodes of interest 

(NODiNT) whose voltages are to be printed and 
plotted against time. The maximum value allowed 
for NODINT is 10. 

FORMAT less. 

6) (IVNCDE (I), I = 1, NCDINT) - This card contains the nodes of 

interest as a request list (IVNCDE), IVNCDE cannot 
include any node that belongs to a subckt. 

FCRMAT less, 

7) NTPIVR,' i^Pl® - This card contains the variable NTPl'iiR, which 

can be either a *Y’ or an *N*, and the total no. 
of elements whose powers are requested to be 
printed end plotted (fCPWR). 

If NTPWR = ’Y' the total power dissipation 
of the circuit is computed and printed/plotted. 

If NTPiiVR = ’N' total power dissipation is not 
calculated. Maximum value of NEPVS/R, allowed is 
10 . 


FCBMAT (3X, Al, IX, 13) 
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8) (IPELM(I,J), 1—1,2), J=l, NEP'//R) — This card contains the 

name (IPELM) (l,J)) and no. (IPELM(2,J)) of the 
elements v^ose powers are to be printed/plotted. 
The request list must not contain any element 
included in a subckt. Skip this card, if mp'm 
in card 7, is zero. 

F0RMAT(3X, 10(A5, IX, 13, IX) 

9) NNXT - This card contains the total no. of nodes (NNXT) 

whose voltages are to be printed in every simu- 
lation step in the file (XJT.NXT. 

FORMAT less 

10) (r^EXT(l), I=1,I^OT) - This card contains the nos. of the 

nodes whose voltages are to be printed in CUT.NXT. 
There must not be > 15 such nodes in the request 
list, 

FCRIi/lAT less. 

Skip this card if NNXT, in previous card is 
equal to zero. 

11) INPRE - This card contains the total no. (INPRE) of 

EXDEF (table 3.2) inputs. 

FOi'vN'A.l less. 

12) (NPREV(I), 1=1, INPRE) - This card contains the nos. of 

those nodes, where EXDEF inputs are applied 


F OFi MAT less. 
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13) PMAX, PIviAXT — This card contains the maximum limits fox 

the power scales in the plots for i) element power 
dissipation (PMAX) and ii) total power dissipation 
(PMAXT). 

FORM! less 

This marks the end of data inputs through files. 

B,3 Interaction though TTY: 

Vi/hen the simulator is run with above data files the 
user may be asked to provide some more data interactively 
through TTY as follows; 

1) If after circuit initialisation, some nodes are found 
floating, (such as output, terminal of an ’off’ pass transistor), 
the user is asked to assign a logic state >0’ or ’I’-to these 
nodes to his liking, 

2) Before voltage/current initialisation, the program 
accepts few NMOS device parameters from the user, through TTY. 
Any parameter typed as ‘O’, is assigned its default value. 

P arameter D efaul t value 


VDD (power supply voltage) 

MU (surface mobility of electrons) 
EPS (dielectric constant of Si02 ) 
TOX (oxide-thickness) 


5V 

800 cm^/V-sec. 
0.345x10'”^ pF/cm 
1.0 X lO”^ cm 
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D efault value 


RATIOD (W/L ratio of driver transistor) 2 

RATIOL (V!//L ratio of . load trans.) 0.5 

RATIOP (vV/L ratio of pass trans.) 1 

RATN2 (W/L ratio of load trans, in a 

two-input' NAiiD gate) ,25 

RATiNl4 (V«/L ratio of load trans, in a 

four-input N/\ND gate) ,125 

CGS (Gate to source capacitance of an 

NMOS transistor) .01 pF 

CGD (Gate to drain capacitance of an 

NMOS trans.) .01 pF 


3) Specifications for the INSEN, CLOCK and PHAS1/PHAS2 
inputs (Table B.2) are accepted through TTY, from the user, 
just before the analysis of the initialised circuit. 

First the nonclock inputs are taken up one by one and 
the user is asked to specify the following for each input; 

a) The input-type (Table B.2) 

b) Timings for each zone (rise/fall edge, high/low 
period) of the input waveform (if the input type 4 - 

’COIST*). 

All timings are in the form NXDELTA, where N is a positive 
integer and DELTA is the simulation time step. To Specify a 
zone timing, the user simply types in the value of N. For a 
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For a rising/falling edge the rise/fall time can not be 
specified less than 4xDELTA. 

After all the inputs have been specified, the 

simulator accepts the timings of the clock inputs, namely 
’CLOCK’ and PHAS1/PHAS2, in the same way as timings for non 
’CONoT’ INGEN inputs are accepted. 

Finally the user is asked to specify the simulation 
period, which is also of the form NXDELTA, 

At the end of the simulation period, the user is asked 
to guide the future course of action of the simulation. At 
this point, 

i) Typing an 'S', stops the simulator 

ii) Typing an ’ R’ , restarts simulation 

iii) Typing a ’ C , c ontinues simulation 

On RESTART’ ing, the circuit is reinitialised and the user has 
to respecify all data, which are accepted by the simulator, 
through TTY. 

On COi'JTINUE’ ing, circuit simulation is carried on, with 
the currently existing state of the circuit as initial state. 
Before starting analysis, the user is asked to type-in the new 
simulation period, ’NPER’ and waveform specifications of inputs, 
if any. The simulator asks the user whether there is any 
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input to be specified in the new period. If there is any 
such, the user types-in the 'input no.* and subsequently gives 
the Vi/aveform specifications, for that input, in the same way 
as told in 3(b). Whereas, if there is no (more) input to 
specify, the user types-in a O and the simulator starts 
analysis over the new simulation period, forth with. Note, 
only non clock IMGEN Inputs can be specified, here. Also, no 
INGEN input can be specified, which has not reached the final 
steady state, at the end of the previous simulation period. 

B.4 GENERAL PROCEDURE: 

V<ith the input formatting as elaborated in above 
sections, the general procedure to be followed, is as 
described here. 

i) The test circuit should be strictly an MMOS digital 
circuit, which can be described with the MOSIMR elements 
(Table B.l), There should be no wired connection of gate 
outputs, except those involving pass transistors. 

ii) Visualise and draw the circuit diagram in terms of 
elements, supported by MOSli^tfl. Interconnections of these 
elements are the circuit nodes and interconnections of devices 
inside an element are the internal element nodes. 
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Number the circuit nodes, sequentially, starting 
from -1. Number elements of each type separately, starting 
from 1, 

iii) If che test circuit is a repeatitive structure, 
where particular blocks can be identified to be connected 
regularly; each individual block can be described ae a 
subckt. Each subckt. is an independent entity where the 
element and node numbering is l ocal , done in the same way as 
in (ii), irrespective of the main circuit numbering. However, 
the nodes in a subckt, are numbered sequentially in such a 
way that internal nodes are numbered first and the external 
nodes (having connections with the environment) are numbered 
later, with higher integers than the internal nodes, 

A subckt, may be located in different parts of the main 
circuit. It is CALL'ed the no. of times it occurs in the main 
circuit, supplying relevant external nodes (relevant Vi/ich the 
main circuit, environment, where it is to be expanded) in each 
CALL. A subckt, cannot include another subckt, 

iv) For proper linking of subckts, with main ckt. 
environment care must be taken in preparing the data in cards 
7 and 9( c) in Sec. B.2. 

Mark the variable LIN (card 7) v/ith an 'L', in a main 
ckt. elem. data for the corresponding input/output node ’N’ 
in the following cases - 
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a) If N is an output node of the element. 

b) If N is an input node of the element and the element 
is the leading member of the NB chain, at N, 

Include a subckt. element in EXTERN (card 9(c)), 
associated with a dummy external node N, only in the following 
cases - 

a) N is an output node of the subckt. element. 

b) N is an input node of the subckt. element and the 
element is the leading member of the NB chain, at N. 

v) EPSLN is the lower limit, set for comparison of a 
node voltage in consecutive tin^ steps, A change in node 
voltage is recognised in the interval t^^ to t 2 , if 

i^n^^l) - ^ 

else V is considered to be constant in that interval, 
n 

This EPSLN and the simulation time step DELTA, are tl^ 
two parameters which can be varied to trade accuracy and 
numerical stability of analysis with speed. 

As EPSLN/^ SPEED and stability and accuracy 

As DELTA /^S>PEED^/^ and stability and accuracy\^ 
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Default values of the parameters (appropriate for node 
capacitances of the order of *01 pF), EPSLN = .00001 V and 
delta = .1 nsec, have been specified in the DATA statement, 
at the beginning of the main pgm. 

If any change in these parameter values are to be 
made the program MOSl^fl is copied and the DATA statement 
in the copy is m-odified, 

vi) Inputs can be applied only at main circuit nodes. 
Nodes and elements to be monitored, must not belong tp a 


subcircuit 
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Table B.l; 


Name of Implication Element specification format 

elem. 



INVER Inverter nt NX 

o ^ 


’ liNlVER* m 


FCRAmT(3X, A5, IX, 13} 


N1 N2 Name of No. of 
NB elem. NB 
at N1 elem. 

at N1 


Name of No. of 
of succ. succ. elem 
elem, at at N2 
N2 


FCBMT(3X,2(I3,IX),2(A5,1X,I3,1X) 

I PAST pass 

transistor o "j o/ntS 

» ! i £ 



‘I PAST* NE 


F Oll'/AT ( 3 X , A5 , IX , 1 3 ) 

N1 N2 N3 Name of No. of Name of No. of 

NB elem. NB elem. NB elem. NB elem. 
at N1 at N1 at N2 at N2 


Name of 
succ .elem. 
at N3 


No. of succ. 
elem. at 
N3 


FCfiivlAT(3X,3(l3,lX),3(A5,lX,l3,IX] 
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Name of ~ 

elem. Implication 

Ype_ 

NAND2 Two-input 
gate 


‘NAND2’ m 
FCBMAT(3X,A5,1X,I3) 

Nl N2 N3 Name of No. of Name of NB 
NB elem, NB elem, elem. 

at Nl at Nl at N2 

No, of Name of No, of 
NB elem, NB succ. succ. elem, 

at N2 elem, at N3 

at N3 

FCaMAT(3X,3(l3,lX),3(A5,lX,I3,lX)) 


N0RG2 Two-input 
NOR gate 


•N0RG2‘ 

FCEiUT (3X, A5, IX, 13) 

Nl N2 N3 Name of No. of Name of 
elem. NB elem, NB elem, 
at Ml at Nl at N2 

No. of Name of No. of succ. 
NB elem, succ, elem, at 

at N2 elem. N3 

at N3 

FCRHAT(3X,3(l3,lX),3(rt5,lX,l3,lX)) 



Element specification format 



Table B.l contd... 
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Hame of 

elem. Implication 

type 


NAND4 


Four- input 
NAND gate 


Element specification format 






NE 



;3x 

, A5, 

IX, 

13) 

Name of 

No 

. of 

NB 

elem. 

NB 

elem 

at 

N1 

at 

N1 

Name of 

No 

. of 

NB 

elem. 

NB 

elem 

at 

N4 

at 

N4 

Name of 


No. ( 

SUi 

cc. elem. 

elem 

at 

N5 




NORG4 


Four-input 
NOR gate 


F0RMAT(3X,5(I3,1X),5(A5,1X,I3,1X)) 

jN 3 


NE 









»N0RG4’ me 
FOR^ iAT(3X,A5,lX,l3) 
N1 N2 N3 N4 N5 


Name of 

No, 

. pf ... 

NB elem. 

NB 

elem. 

at N1 

at 

N1 

Name of 

No 

. of 

NB elem. 

NB 

e lem 

at N4 

at 

N4 

Name of 


No. of succ. 

succ. elem. 

elem. at 

at N5 


M5 


Table B. 1 contd.. 
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Name of 

elem. Implication Element specification format 

type 

NI3UF Nohinverting 
super buffer 

’NIBUF» NE 

FORMAT (3X, A5, IX, 13) 

N1 N2 Name of No, of Name of No. of 

NB elem. NB elem. succ. succ. 

at N1 at Ml elem. elem, 

at N2 at N2 

F0RMAT(3X,2(I3,1X),2(A5,1X,I3),1X)) 


INBUF Inverting 

super buffer 


' INBUF ' NE 

FORMAT (3X, A5, IX, 13) 

N1 N2 Name of No. ofiMB Name of 

NB elem. elem. succ. elem. 

at N1 at 'Nl at N2 

No, of succ, 
elem. at N2 

FORMAT (3X,2(I3,lX),2(A5,lX,I3,iX)) 




Table B.l contd.,. 
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Name of 
elem, 

LATCH 


‘LATCH’ ^E 

FCBMAT (3X, A5, IX, l3) 

N1 N2 N3 N4 Name of No. of Name of !© 

NB elem. NB eleia. elem* 
at N1 at N1 at N2 

• No. of Name of No. of 
MB elem. succ. Succ. 
at N2 elem. at elem.at 
N3 N3 

Name of No. of 

succ, elem. succ, elem. 
a t N4 a t N4 

F(HMAT(3X,4(I3,1X) ,4( {A5,lX,l3,iX)) 

NOELM No element ‘NC3ELM’ OCX) 

FORMAT (3X, A5, IX, I3) 

■ If a subckt, elem. has a main circuit 
element, as its NB/succ. or a main ckt. 
elem. has a subckt, elem. as fts NB/succ., 
the corresponding linking is flone during 
subckt. expansion, v/hen the subckt is 
CALL'ed. But in the circuit description 
(the subckt. and main ckt. being independent) 
the corresponding NB/succ. field of data 
(cards 5 and 8) is filled with ‘NOELM*. 




N.B. NE = element no 
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Table B. 2: Input Table 

General format for any : 'INPUT* NI 

kind of input FCRMATCSX, Ab, IX, 13) 

NA Input kind Name of No. of succ, 

succ. elem. 
elem. 


FCaMAT (3X,I3,1X,2(A5,1X)J3) 


I nput kind Implication 


specifications 


EXDEF Externally defined The voltages at E>DEF input nodes 

input are specified in INP.PRE in 

proper order, at each simulation 
time step, in the format 

FCEMAT( 3X, 15(F7. 3, IX) 

over the entire simulation period, 

INGEN Internally Timings such as 

generated input 

ion after which input 
° * commences to change 

1) REDGE(lislng edge) Rise/fall time 

ii) FEDGE(faliing iii) High/low period are 

fed-in interactively 

iii) PRJLS( positive through TTY 

pulse) 

iv) NPULS( negative 
pulse) 

v) CONST( constant) 

CLOCK single phase clock 

PHASl/ Two phase clock 
PHAS2 


N,,B.; NI = Input no. 

NA = Node where input is applied. 
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